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ABSTRACT 

A comprehensive experimental and analytical investigation has been 
undertaken to determine the hydroelastic response of elastically-mounted 
smooth and rough cylinders in harmonically-oscillating flow. 

The forces acting on the cylinder and the cylinder response have 
been measured together with the characteristics of the flow and analyzed 
through the use of the Fourier analysis, spectral analysis, Duhamel's 
integral, and the wake-oscillator model. In addition the in-line force 
has been analyzed through the use of Morison's equation. 

The results have shown that hydroelastic oscillations occur ina 
narrow range of the reduced velocity and give rise to significant 
amplification of the exciting force. 

The use of the Fourier and spectral analysis together with the 
governing equations of motion was shown to predict the hydroelastic 
response of the cylinders in the synchronous region. 

The investigation shed considerable light on the interaction between 
the harmonic fluid motion and the hydroelastic response of cylinders and 
provided definitive information for the prediction of the force and 


response characteristics of the type of oscillations considered herein. 
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IT. GENERAL SURVEY 


A. INTRODUCTION 

Experience has shown that elastic structures near linear-resonance 
conditions can develop flow-induced oscillations by extracting energy 
from the flow about them. The oscillations of the structure modify the 
flow and give rise to nonlinear interaction. This is in addition to 
any non-linearity which can arise from the restoring force (variable 
support stiffness) and/or from response-dependent structural damping. 
The understanding of these nonlinear interactions is of paramount importance. 

Four general classes cover a wide range of the fluid-structure 
interaction phenomenon: vortex-induced oscillation, flutter, galloping 
and buffeting. This survey will deal with only the first of these 
phenomena. 

The problems associated with vortex-induced oscillations and their 
important consequences have proved to be exceedingly difficult both 
theoretically and experimentally. The difficulty is partly due to our 
incomplete knowledge of the flow field about a stationary body but it is 
also due to fundamental problems associated with the coupling of the 
oscillations and the flow field. It would indeed be gratifying to note 
that our understanding of the separated flow about a bluff body, particu- 
larly in the region enclosing the body and its near wake, is complete 
through laboratory and numerical experiments. This is, regrettably, not 
yet the case, but there have been nevertheless significant advances since 
the time attention has shifted from idealized vortex streets to the 


details of separation, three-dimensional flow, and the role played by the 
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afterbody on the formation, growth, and motion of vortices. Many of the 
key developments of the past decade have been incorporated into the 
experimental investigation of flow-induced oscillations and into the 
evolution of conceptually plausible idealized models. 

In the material to follow emphasis is placed on the understanding of 
the transverse oscillations of cylinders. No attempt has been made to 
be encyclopedic in coverage as there are recent texts and survey articles, 
for example, Blevins [1], Parkinson [2], King [3], and McCroskey [4], 


which cover the field in some detail. 


B. VORTEX SHEDDING FROM A STATIONARY BLUFF BODY 

There are substantial gaps in the understanding of the fundamental] 
mechanisms at work in the forcing of the body by the forces acting on 
its afterbody and in the prediction of the lock-in range in which resonant 
oscillations occur. To a degree, the difficulties in obtaining a solution 
or in devising a model stem from an inadequate knowledge of the precise 
mechanism of flow separation and its cOnsequences in steady flow past a 
stationary body (see, e.g., Morkovin [5], Marris [6], Mair and Maull [7], 
and Berger and Wille [8] for comprehensive reviews and references) . 

Observations as well as numerical experiments have shown that the 
wake of a bluff body is comprised of an alternating vortex street. This 
phenomenon is intrinsic to the flow itself and is a consequence of the 
interaction between the shear layers, base pressure, diffusion and dissipa- 
tion of vorticity, and far wake, i.e., a consequence of the particular 
vorticity distribution throughout the field. 

As the vortices grow symmetrically, say, in an impulsively-started flow 


about a cylinder, the shear layer joining the separation point to one of 
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the vortices begins to develop instabilities (presumably Tollmien- 
Schlichting instabilities [9]) and is drawn across the wake in response 

to the base pressure reduced by the action of the vortex growing across the 
wake (see Fig. 1) [10]. This nearly corresponds to a time when the sheet 
drawn in has least circulation or is most permeable. The stretching, 
diffusion, and dissipation of vorticity break up the deforming turbulent 
sheet and thereby the further supply of circulation to the vortex whose 
rate of growth has already been reduced to its minimum. This corresponds 
to the shedding of the first vortex. The vortex across the wake still 
continues to grow (but at a decreasing rate) and entrains part of the 
Oppositely-signed vorticity left in the wake by the cut sheet and the 
irrotational fluid drawn from outside through the opening created by the 
shedding of the first vortex. The shedding process for the second vortex 
does not commence until the circulation in its feeding sheet decreases to 
its minimum, making the sheet most susceptible to rapid diffusion. Simul ta- 
neously, the sheet deforms, diffuses, and is drawn across the wake by the 
action of the base pressure and the vortex growing on the other side of the 


The above mechanism, 


wake. Then the shedding cycle repeats itself. 
exhibited by numerical experiments [10] is quite similar, but not identical, 
to that suggested by Gerrard [ll]. 

The separation point at the side of the cut sheet just passes through 


its minimum angle. The sheet at the opposite side of the wake bears nearly 


the maximum circulation, relative to any other time, and its separation-point 


Ve or a circular cylinder there is a 60 degree phase difference between 
the lift maximum and the sheet cutting time, i.e., a vortex is shed when 
the lift force decreases to 50 percent of the absolute value of its maximum. 
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Fig. la. Initiation of vortex shedding [10] 
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angle just passes through a maximum. The motion of the stagnation point 
is such that it is 180 degrees out of phase with the separation points, 
i.e., the separation point is below the x-axis (flow from left to right) 
when a vortex is shed from the bottom of the cylinder. 

The shedding of the vortices gives rise to an oscillating side thrust, 
upon a cylinder of suitable afterbody, in a direction away from the last 
detached vortex. The relationship between the vortex-shedding frequency, 
cylinder diameter, and the velocity of the ambient flow has been discovered 
by Strouhal [12] almost exactly 100 years ago in connection with his work 
On a Special method of creation of sound. His data showed that the product 
fD/V is nearly constant for a wide range of values of D and V. This pos- 
tulated constant is often called the Strouhal number and Strouhal's own 
data suggested a value of about 0.185. The characterization of the vortex 
shedding process by a simple frequency is a practical simplification. As 
first noted by Lord Rayleigh [13], the Strouhal number is a function of the 
Reynolds number for a given body? and f should be written as f = S(Re)V/D, 
(see Fig. 2). 

At relatively low subcritical Reynolds numbers, the energy containing 
frequencies are confined to a narrow band, and the Strouhal number is about 
0.20 for smooth cylinders [8]. For Reynolds numbers larger than about 
20,000, only an average Strouhal number may be defined. In the range 


0.5x10° 


< Re < 2.5x10° no periodicity is observed in either the unsteady 


4s, = 0.21(1 - 20/Re) for 40 < Re < 1,000 (Goldstein [14]); 
S. = 0:212(1 - 12.7/Re) in the range 400 < Re < 100,000 according to a 
cOrrelation obtained by Roshko [15]. In the range 400 < Re < 1,200 the 
results of Roshko's correlation are slightly below the best-fit line to 
the existing data [8]. 
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lift or the unsteady pressures. In the range 2.5x10° < Re < 6.5x10° 


a small peak is observed at 5, = 0.2 in the unsteady lift and in unsteady 
pressures only at the front part of the cylinder. No peak is found behind 
the separation point. For Re > 6 .5x10°, a definite peak is found in both 
the unsteady lift and the unsteady pressures at So = 0.3. The said 
periodicity occurs both at the front and rear part of the cylinder [16]. 
The foregoing categorization is somewhat arbitrary. Jones at al., 

[17] classified the unsteady lift force in three regimes: wide-band 


6 Re < 3.5x10°; narrow-band random for 3.5x10° <aRe < 6x10": 


random for 1.1x10 
and quasi-periodic for 6x10° < Re < 18.7x10°. 

The spectral content of the exciting forces is particularly important 
for bodies which may undergo in-line and/or transverse oscillations since 
the vortex shedding frequency locks on to the frequency of the transverse 
oscillations of the cylinder when there is a transverse force and when the 
vortex shedding frequency is in the neighborhood of the natural frequency 
of the cylinder. When the flow about a freely oscillating cylinder is in 
the critical regime, the sharp rise in the Strouhal number does not occur 
and remains at a value nearly equal to that found at subcritical Reynolds 
numbers [18]. 

Various attempts have been made to devise a universal Strouhal number 
which would remain constant for differently shaped two-dimensional and 
axisymmetric bodies (see Table-I). In Table-I h' represents the wake 
width as determined from Roshko's notched hodograph theory [19]; Qs 
the longitudinal spacing of the vortices as determined by Kronauer's 
minimum drag criteria as reported by Bearman [20]; Qs the wake width at 
a distance of the wake establishment region; and KY = (1 - A 6 i 


where Cob is the base pressure coefficient. 
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Table I Universal Strouhal numbers 


So value Body shape Source 
Fh /VK 0.163 two-dimensional Roshko [19] 
Fe y/V Oris two-dimensional Bearman [20] 
By! “Ky 0.190 axisymmetric Calvert [21] 
Fey! VK, 0.163 two-dimensiona | Simmons [22] 


The foregoing definitions suffer from the obvious drawback that they 
require the solution of the wake formation problem first, or the 
measurement of one or more flow characteristics. 

Theoretical or semi-empirical predictions of the wake characteristics 
have been attempted by various researchers. Birkhoff [23] demonstrated 
that the longitudinal spacing is trivially invariant since the longi- 
tudinal velocity of vortices is uniformly bounded. He has shown further 
that in an inviscid fluid the lateral spacing also remains constant at a 
value equal to the initial spacing of the shear layers. These led Birkhoff 
to the conclusion that h/k,, is determined by its initial value and not 
by von Karman's stability criteria of h/d,, = 0.281. The use of a wake 
oscillator model ("the wake swings from side to side, somewhat like the 
tail of a swimming fish" [(23]) and some experimentally-determined values 
led Birkhoff to So = 0.2. One must keep in mind that vorticity in real 
vortices is not concentrated in points, the vortices are non-circular and 
distort and rotate as they move downstream [24], vorticity diffuses and 
is swept across the wake [25], and dissipates by turbulence [8]. Vortices 
are subjected to strain fields imposed by nearby vortices. The resulting 


patterns are ever changing vortex shapes encompassing elliptic and 
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pearlike geometries. The complexity of the interaction between strained 
distorting vortices is further exacerbated by the addition of turbulence 
to the wake, as this is likely to produce a more diffusive vorticity 
distribution and thus an additional shear field [24]. Thus, the near 
constancy of wake characteristics in the range of Reynolds numbers where 

a vortex street might be observed is primarily due to the slow variation 
of dis hy» and So with respect to the strength and deformation of vortices. 
Nevertheless, theoretical idealizations of the wake give some clues about 
the asymptotic behavior of the wake. In particular, one obtains from 


Fey = \V - Us where ee is the velocity relative to the body, and 


from the use of von Karman's stability criteria that 


1 P/VO 


Sa. = [= Ol 
Oo RAD 2v2(2, /D)* 





where [T represents the strength of a vortex. Assuming that &/D remains 
nearly constant (2,/D = 5), one concludes that the smaller the strength 

of the vortices, the larger is the Strouhal number (smaller drag and 
narrower wake) and vice versa. The vortices which are subjected to 
smaller dissipation linger a little while longer in the near wake relative 
to those which are dissipated more, i.e., strong vortices get stronger. 
However, the dependence of 36 on F/VD-is rather weak. This explains in 
part the reason for the success of the discrete vortex models in predicting 
the Strouhal number fairly accurately [26,27] in spite of the fact that 
the calculated vortex strengths are about 35 percent larger than those 
estimated experimentally. However weak, the dependence of So on T/VD 
becomes a primary factor in the locking-on of the vortex shedding to the 


natural frequency of the body. The vortices which are stronger continue 
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to be fed by their shear layers a longer time period, thus further reducing 
the Strouhal number. In this process the mobility of the separation 
points is important but not necessary. 

Sacksteder [28] pursued a theoretical approach to determine the Strouhal 
number at large Reynolds numbers by perturbing the d'Alembert flow and 
obtained So = 0.2028. His analysis does not consider the wake formation 
and thus it is not expected that flows around objects that induce large 
wakes could be treated with a simple perturbation of the unseparated flow. 

The constancy of the Strouhal number over a broad range of Reynolds 
numbers does not imply that the base pressure remains constant and that 
a simple two-dimensional vortex emanates from a separation line. In 
reality, there is not only a phase shift between various sections of the 
vortex, separated by a correlation length,> but also variations in both 
the intensity and the frequency of vortex segments. 

The variation of the base pressure with Reynolds number, in the range 
where the Strouhal number remains practically constant, may be related to 
the variation of the mean vorticity flux or to the variation of the 
correlation length with the Reynolds number, turbulence, length-to-diameter 
ratio, and surface roughness [29-33]. Table-2, as compiled by King [3], 
gives an approximate idea about the typical values of the correlation 
length. The net effect of the spanwise variations of the vortex tube is 
that the transverse force (lift) coefficient obtained from a pressure 
integration is not necessarily identical with that obtained from the direct 


measurements of the lift force. Partial spanwise correlation leads to 


3The equivalent length over which the velocity fluctuations in the 
wake may be described as perfectly correlated. 
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Table 2 Correlation lengths 


Reynolds No. range Correlation length Source 
40 < Re < 150 15D - 20D [34] 
150 < Re < 10° 2D - 30 [34] 
10° < Re < 4.5x104 30 - 6D [32] 
Re > 10° 0.05D [34] 
Re = 2x10° 1.56D [29] 


variations in both the frequency and the amplitude of the lift force, the 
variation of the latter being more pronounced than that of the former. 
The lack of correlation exists not only spanwise but also chordwise 
[35-38] and the chordwise correlation is related to the spanwise correla- 
tion. Comparison between various results suggests that with increasing 


4 to 10°, the chordwise correlation for 


Reynolds number over the range 10 
Square and circular cross-section cylinders is maintained or improved. This 
leads to an increase in fluctuating lift coefficient. The reasons for 
these variations are not quite clear. The end effects [29, 39-41], wall 
boundary layers, freestream turbulence [42-44], and non-uniformity of 
the flow are mentioned often as possible reasons. The complexity of the 
three-dimensional nature of the flow about a cylinder is clearly demon- 
Strated with measurements by Tournier and Py [45]. 

It would not be correct to assume that the mobility of the separation 
points is primarily responsible for the imperfect coherence. Even bodies 
such as 90-degree wedges, square cylinders, with fixed separation lines, 


do not exhibit perfect correlation. However, the variation of the base- 


pressure coefficient for bodies with mobile separation lines is greater 
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than that for bodies with fixed separation lines [46]. In fact, for 
flow over symmetrical wedges the base pressure appears to be insensitive 
to the nose angle [46]. 

As noted earlier, the strength of the vortices plays an important 
role, particularly in the near wake. Laboratory and numerical experiments 
have shown that [/, 19, 47-50] the net circulation of a rolled-up vortex 
of the street is 40 to 60 percent smaller than that generated in the boun- 
dary layer during a shedding cycle. Prandtl] determined that the initial 
vorticity decreases to about half where the first vortex centers appear. 
Vorticity is ultimately dissipated by viscosity to which it owes its 
generation. Nevertheless, one may think of loss of circulation through 
cancellation of oppositely-signed vorticity. Primarily, there are three 
mechanisms whereby oppositely-signed vorticity are brought close together: 
vorticity generated on the forebody is carried by the shear layers near 
that generated on the afterbody; vorticity of the deformed and cut sheet 
is carried across the near wake by the entrainment of the irrotational 
fluid; and, finally, vorticity is swept across the entire wake [25]. The 
percentages quoted in the literature [8] for the total loss of vorticity 
often imply that the vortices, once having acquired a certain circulation, 
retain that circulation throughout the rest of their motion. The fact 
that circulation decreases continuously with time or distance is demonstrated 
clearly by the experiments of Schmidt and Tilmann [50] and Bloor and 
Gerrard [48]. The amount of vorticity generated in the boundary layers 
and the amount dissipated are of prime importance not only for the flow 
about stationary bluff bodies but also for those undergoing resonant 
oscillations. In fact, the entire bluff body problem may be reduced to 


the determination of the vorticity distribution throughout the flow field. 
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This is not yet possible for Re larger than about 100. The determination 
of the vortex strengths is difficult and sensitive to the theoretical 
and experimental means employed. 

It is evident from the foregoing and from a more detailed perusal of 
the references cited that the description of the near wake of a bluff 
body is in a primitive state. Much of what is known about the consequences 
of separation has come from laboratory experiments. It has not yet been 
possible to develop a numerical model with which experiments may be con- 
ducted to explain the observed or inferred relationships between various 
parameters and to guide and complement the laboratory experiments. The 
principal difficulties are as follows: (i) Separation points. They repre- 
sent a mobile boundary between two regions of vastly different scales. 
This in turn leads to complex physical nonuniformities in relatively 
narrow regions which cannot be handled within the framework of the boundary- 
layer theory [51]. Finite difference and Marker and Cell (MAC) techniques 
require in such regions very smal] computational times. The discretization 
of the continuous process of vorticity generation by line vortices in the 
vicinity of a mobile or fixed singular point (discrete vortex model ) 
strongly affects the existing non-uniformities and promotes earlier 
Separation. Attempts to preserve the prevailing flow conditions, say by 
limiting the influence of the nascent vortices, while satisfying a rela- 
tively simple separation criterion lead to hydrodynamical inconsistencies 
and non-disposable parameters. (ii) Reynolds number. Finite difference 
schemes for bluff-body flows are limited to relatively small Reynolds 
numbers whether the scale of flow is assumed to be governed by a constant 
Viscosity or by a constant eddy viscosity. The large recirculation region 


of the flow is often comprised of turbulent vortices even when the boundary 
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layer is laminar. The transition to turbulence moves upstream in the 


shear layers as the Reynolds number is increased from about 10° to 


0. At Re = 5x10", it reaches the shoulder of the cylinder [52]. 

It does not move appreciably further upstream before the critical Reynolds 
number is reached. Thus, the distribution, turbulent diffusion, and 

decay Of vorticity and the interaction between the wake and the boundary 
layers cannot be subjected to numerical simulation without recourse to 
some heuristic turbulence models and inspired foresight. 

The representation of the wake by clouds of point vortices or dis- 
cretized spiralling sheets (see, e.g., [10, 26, 49, 53]) is not immune to 
scaling problems. In fact, not a particular Reynolds number but only a 
particular flow regime may be specified, depending on the separation cri- 
teria used. (iii) Three-dimensionality. As noted earlier, even a uniform 
flow about a stationary cylinder exhibits chordwise and spanwise variations. 
These three-dimensional effects may play a major role in the stretching of 
vortex filaments and in the redistribution of vorticity in all directions. 
The numerical models are not in a position to account for such complex 
effects. Qne may hope to assess the effects of three dimensionality by 
means of two-dimensional numerical experiments. 

It is against this background that the nonlinear fluid-structure 
interaction projects an order of magnitude more complex problems. The 
picture emerging from the numerous experiments and mathematical models 
is that there are phenomena to be explained, parameters to be quantified, 


and conjectures to be justified. 


C. LOCK-IN OR SYNCHRONIZATION IN STEADY FLOW 
Numerous experiments have shown that when the vortex-shedding frequency 


brackets the natural frequency of an elastic or elastically-mounted rigid 
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cylinder with a suitable afterbody (capable of giving rise to a trans- 
verse force), the cylinder takes control of the shedding in apparent 
violation of the Strouhal relationship. Then the frequencies of vortex 
Shedding and the body oscillation collapse into a single frequency close 
to the natural frequency of the body (see Fig. 3). This phenomenon is 
known by various names: lock-in, locking-on, synchronization, hydro- 
elastic or fluidelastic oscillations, wake capture, self-controlled or 
self-excited oscillations, etc.* 

The facts which have emerged from two decades of work on vortex- 
induced oscillations may be summarized as follows: 

1. When a body is close to its linear resonance conditions it can 
undergo sustained oscillations at a frequency close to its natural 
frequency [55]. That a cylinder should be excited at its natural fre- 
quency when the frequency of the exciting force is equal to its natural 
frequency is not surprising. But that the phenomenon encompasses a range 
of +25 to 30 percent of the natural vortex-shedding frequency and that the 
vibration and vortex shedding frequencies lock together and control the 


Shedding process are surprising. 


ui vortex-excited oscillation is actually a forced one having a 
self-excited character also to some degree due to lift force amplifica- 
tion through nonlinear interactions. By definition, sel f-excited 
oscillation is one where the alternating force that sustains the motion 
is created or controlled by the motion itself; when the motion stops, 
the alternating force disappears. For the phenomenon discussed herein 
the alternating lift force does not disappear when there is no oscilla- 
tion. In fact, when there is no alternating force (a flat plate normal 
to the flow) there is no lock-in. Apparently, the often-used definition 
of 'self-excited oscillation’ is a misnomer. 
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Fig. 3 Response characteristics of a freely-oscillating 
cylinder, from Feng [54], (see also [2]) 
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2. The interaction between the oscillation of the body and the 
action of the fluid is nonlinear (e.g., [56-57]). 

3. The vortex shedding does not necessarily result in an alternating 
transverse force. There must be a suitable afterbody, and hence an 
alternating lift force, and the shedding frequency should coincide with 
or be near the natural frequency of the body, or one of its harmonics, 
for sustained oscillations to occur. Thus, the magnitude and occurrence 
of sustained oscillations strongly depend on the lift coefficient of the 
Stationary body. 

4. The parameters ve = V/FD and So = mc/oD° [58] or RO me/pLD“(C) 
[60] or A, = (2nz)(nS,)°/p, are of major importance in determining the 
amplitude of oscillations and the range of synchronization for a given 
body [59-61]. 

5. The velocity range over which synchronization occurs increases 
with oscillation amplitude [62]. 

6. The excitation range of cross-flow oscillations in air extends over 
4.75 < ve < 8 and maximum amplitudes occur in the range 5.5 < ve <5) 

[3]. In water, the excitation range of the transverse oscillations can 
be increased to 4.5 < ve < 10 with maximum amplitude falling within the 
range of 6.5 < _ <6. 

7. For a circular cylinder with large L/D, synchronization begins 
when lig ie and ends at about Ue ~ 1.4. The maximum amplitude occurs 
near the middle of this range. At the end of the lock-in range, vortex- 
Shedding frequency jumps to that governed by the Strouhal relationship, 
but the cylinder continues to oscillate at f ~ f,. This is true at 
both ends of the lock-in range [54] and shows that the response is not a 


Simple forced vibration at the exciting natural Strouhal frequency. 
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8. For a cylinder with fixed separation points (e.g., a D-shaped 
cylinder) synchronization extends over a range of 0.96 to 1.1 of the 
resonant velocity Sele from about 0.85 to 1.1) [24]. According to 
Feng [54], the said range extends over a range of 5 < Va < 7. Gowda [39] 
found that the lock-in begins already when ii is about 80 percent of ae 
and ends when i =f. The maximum amplitude occurs near the end of the 
synchronization. 

9. In synchronization, the slantwise vortex shedding is replaced by 
parallel vortex shedding [44, 62]. 

10. The correlation length increases rapidly with amplitude [57, 62-65]. 
The increase of the correlation length in smooth flow is much larger than 
in turbulent flow [66]. In smooth flow the correlation length is estimated 
by numerical extrapolation to increase from about 3.5D to 40D for 
Re = 1.9x10" In the range” 0.05 < A/D < 0.1. In turbulent flow, it is 
again estimated to vary from about 2.5D to 10D in the same A/D range. 

The rate of increase is steeper than linear but does not show any abrupt 
change which would indicate a sudden development of the lock-in once a 
threshold amplitude is achieved. 

11. For bodies, with mobile or fixed separation points, undergoing 
sustained oscillations, the vortex strength is increased [24, 67]. This 
could be either or both due to increased rate of vorticity flux or due 
to less destruction of vorticity in the near wake. The growing vortex 


On a vibrating cylinder seems to roll up more quickly [24]. 


“The generality of these values is not certain since in Novak and 
Tanaka's [66] experiments the cylinder length (between the end plates) 
was Only 12 diameters. 
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12. The transverse force needed to excite a cylinder to large ampli- 
tude oscillations is far greater than that exerted by vortex shedding, 
i.e., the transverse force is amplified (e.g., [54], [56], [61]). The 
increase of the correlation length is responsible only for about 10 
percent of the lift increase, i.e., it is not the major factor for the 
lock-in phenomenon. 

13. A hysteresis behavior may exist in the amplitude variation and 
frequency capture depending on the approach to the resonance range -- 
whether from a low velocity or from a high velocity [54]. No universal 
behavior is noted and the reasons are not clear. The jump condition 
(double amplitude response) may originate in the fluid system, and there- 
fore in the lift force, and not in the cylinder elastic system [2]. The 
‘combination-oscillation' model [68], (based on the assumption that the 
lift oscillator has two components, one at the cylinder frequency, repre- 
senting lock-in, and the other approximately at the frequency of vortex 
Shedding for the stationary cylinder) gives rise to jumps between two 
Stable oscillation amplitudes but only at the ends of the lock-in range 
and not in the middie of the said range, as found by Feng [54], (see Fig. 3). 
The jump may be the consequence of a variable structural damping (dependent 
on the interaction of the ambient velocity and structure support, vibra- 
tion velocity and/or amplitude, etc.) [69-70] or of a nonlinear spring 
behavior [71]. It is also entirely possible that the jump condition 
originates in both the fluid and the cylinder elastic system, including 
blockage effects [72] and the vibrations of the test apparatus [54]. 

14. The absolute value of the base pressure and the in-line or drag 


force increase with the amplitude of oscillation [3, 56, 61, 72, 73]. For 
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a circular cylinder, as increases as much as 60 percent for A/D = 0.29 
(721.° For a D-shaped cylinder, Davies [24] found an increase of 70 
percent relative to the stationary D-cylinder value. 

15. The phase difference between the displacement and force signals 
increases almost linearly at first and then rapidly [54, 56]. 

16. The mobility of the separation points is important but not 
necessary for synchronization. The separation point on a circular cylinder 
moves over an arc whose magnitude depends on the frequency and amplitude 
of the oscillation. The magnitude of the separation-point travel reaches 
a peak when the cylinder frequency is about 0.9 of the wake frequency 
ma, 75). 

17. The actual instantaneous value of the wake angle is greater than 
that between the ambient flow velocity and the relative fluid velocity 
(74, 75]. 

18. The overall wake width is practically unaffected in spite of the 
lateral movement of the body that generates it. Koopman [62] inferred 
from flow-visualization studies an apparent decrease in lateral vortex 
Spacing in the resonant wake of an oscillating circular cylinder. Davies 
[24] found about a 20 percent increase for a D-shaped cylinder synchronized 
in a flow at much high Reynolds numbers (7x10° < Re < 4x10°, py Gee 2) . 
Griffin [67] found that the lateral spacing of the vortices decreases 
with increasing amplitude and is unaffected by changes in frequency of 


oscillation. 


OF i xing separation at 6. = +70 degrees by a tripping wire on a circular 
cylinder yielded a nearly constant C b at -1.39 for A/D = 0.1 for all 
faues of f /f. from 0.4 to 1.3 (72]P° The reason for this anomalous 
behavior is°no% clear. Confirmatory experiments are needed. 
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19. Longitudinal spacing d does not change with the amplitude of 
motion in the synchronization region but it varies inversely with the 
frequency ratio [77]. For a cylinder undergoing forced transverse 
oscillations with the frequency ee &,, decreases from its value at 
iZ = i. as adh es increases and vice versa. In other words, the wake 
contracts for Waite > | and expands for Tee ssl 

20. The formation length Le exhibits a similar variation with ain 
However, unlike the longitudinal spacing, the formation length decreases 
systematically with increasing amplitude of vibration [76]. For a given 
Fits the vortices in the near wake go through a spacing adjustment from 
an amplitude-dependent spacing to an amplitude-independent longitudinal 
Spacing. 

21. Forced vibrations may differ from the naturally occurring vortex- 
induced vibrations since the latter depend, to some extent, on the history 
of the fluid motion. History effects may be particularly important in 
cases where hysteresis occurs. Driven models overshadow the intricate 
interaction between the evolution of the motion and the variations of 
the pressure distribution. 

22. Frequency demultiplication or frequency division (existence of 
the natural vortex-shedding frequency in the state of synchronization 
in the wake) may occur when the forced cylinder frequency Fe is an 


integral multiple of the shedding frequency’ [S05 7 5 aay | 


For example x = 2b-@sin t is the exact solution of 
* oe 
X - Jaw(1-bx*)x f HOR = -4ab “sin(3ut) 


and is a special case of the phenomenon called ‘frequency demultiplication 
by van der Pol [78]. 
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23. In-line oscillations occur within two adjacent regions [79-86]. 
The first is in the range of 1.25 < ue < 2.5, maximum amplitudes occurring 
at ve = 2.1. The second region extends from Me = 2ato V. = 3.8 with 
maximum amplitudes at Ve = 3.2. The first instability region is accom- 
panied by symmetric vortex shedding (as if the flow started impulsively 
from rest at each cycle) and the second by alternate vortex shedding. 

Several mathematical models have been proposed [87-91], (see also 
[2, 55, 71, and 92) for additional models) in an attempt to simulate 
and/or explain some of the experimental results summarized above. These 
models do not include the analysis of the flow field and the fluid- 
mechanical arguments invoked in their evolution are not altogether con- 
vincing. Thus, their worth should be measured not so much by their 
capacity to obtain functional relations among significant parameters 
that lead to the basic understanding of the phenomenon but rather by 
their ability to produce results which are qualitatively similar to those 
Obtained experimentally. 

The most noteworthy among the oscillator models is the one proposed 
by Hartlen and Currie [87] where a van der Pol type soft nonlinear 
oscillator for the lift force is coupled to the body motion by a linear 
dependence on cylinder velocity. This model is based partly on a sugges- 
tion by Birkhoff and Zarantonello [23], mentioned earlier in connection 
with the estimation of the Strouhal number, and by Bishop and Hassan [56] 
in connection with their now classic experiments with oscillating cylinders 
in uniform flow. The model has its roots in mechanics and electricity 
rather than in the equations of fluid motion. The details of this model 


will be described later. 
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D. LOCK IN OR SYNCHRONIZATION IN OSCILLATING FLOW 

There has been relatively little theoretical and experimental work 
on the complex dynamic response of elastic structures to oscillating 
flows. 

Laird [93] explored, in 1962, the effects of support flexibility by 
oscillating a vertical cylinder through still water. He found that the 
forces acting on a flexibly-supported oscillating cylinder can exceed 
4.5 times the drag force on the cylinder rigidly-supported while moving 
at a uniform velocity equal to the maximum velocity during the oscillation 
and that a cylinder, flexible enough to have transverse oscillations with 
amplitudes more than half the diameter, while performing large amplitude 
oscillations in water, tends to oscillate transversely at the eddy-shedding 
frequency and to vibrate at twice the eddy frequency in the in-line 
direction. Laird's work pointed out the significance of the vibrations 
in harmonic flow and the need for a detailed investigation of the phenomenon. 

Vaicaitis [94] investigated the response of deep-water piles due to 
cross-flow forces generated by wind-induced ocean waves. The resulting 
cross-flow forces were treated as random processes in the time-space 
domain and are assumed to be dependent on fluid velocities and vortex 
Shedding processes. Out of necessity, Vaicaitis had to make a number of 
assumptions regarding the structural motion, the added mass coefficient, 
Strouhal number, etc. Vaicaitis was careful to note that some of his 
conclusions were based on a particular example and that they might not 
be true for piles with different dynamic characteristics. 

Selna and Cho [95] investigated the in-line resonant response of a 
tall structure by assuming constant drag and inertia coefficients in 


calculating the exciting fluid force through Morison's equation [96], 
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modified to take into account the motion of the structure. Their calcu- 
lations have shown that a resonant motion is possible and that the dynamic 
deflection of a flexible structure can considerably areacd its static 
deflection. 

Verly and Every [97] measured wave-induced stress on similar rigid 
and flexible vertical cylinders in a wave channel at relatively low 
Keulegan-Carpenter and Reynolds numbers. Even though they were unable 
to correlate their data with any suitable parameter governing the motion, 
they concluded that the vibration is caused by the cylinder's response to 
eddy shedding and that there is no fluid-structure interaction. They 
found that the vibration occurs if UF D (U_, represents the maximum 
velocity of flow or wave in a cycle; and D, the diameter of the cylinder) 
is greater than about unity for any natural frequency, wave frequency, 
and damping. The vet D parameter never reached high enough values in 
Verly and Every's experiments for the cylinder to undergo synchronized 
oscillations as the present investigation has shown. 

McConnell and Wilson [98] conducted experiments similar to those of 
Laird [93] by sinusoidally driving a circular cylinder in still water. 
These experiments were conducted at rather small Reynolds numbers with 
0.25 in and 0.50 in cylinders. Furthermore, these investigators measured 
only the amplitude of oscillation and calculated the lift coefficient at 
resonance from a slightly-modified version of an equation developed by 
Sarpkaya [60]. They have concluded that et 0 is one of the most impor- 
tant parameters governing the resonant oscillations. 

Sawaragi et al., [99] investigated the in-line and transverse dynamic 


response of a cantilevered circular cylinder, with a concentrated mass 
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at its top, in waves of smal] amplitude. The Reynolds number ranged 

from 1,500 to 6,200 and the rms value of the Keulegan-Carpenter number 

(K = UT/D; calculated over the submerged length of the cylinder through 
the use of rms value of the maximum of the horizontal velocities) ranged 
from 2 to about 20. They have approximated the lift coefficient for a 
rigid cylinder by a Rayleigh distribution and calculated the dynamic 
response of the test cylinder. The results are of limited value since 
the interaction between the synchronization and the force amplification 
was ignored and since no more than three vortices were shed during a half 
cycle. Furthermore, the rapid variation of the Reynolds number, Keulegan- 
Carpenter number, vortex-shedding coherence, and lift force with depth, 
in a range of Keulegan-Carpenter numbers where the phenomenon is least 
understood, obscured the role played by the governing parameters. 

It is clear from the foregoing that there are no well-established 
analytical methods either to analyze the data or to predict the character- 
istics of the flow or of the vibration. Thus, it was decided to undertake 
a comprehensive investigation of the dynamic response of elastically- 
mounted smooth and rough cylinders in harmonic flow particularly at high 


Reynolds numbers and Keulegan-Carpenter numbers. 
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TI. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A. U-SHAPED OSCILLATING-FLOW TUNNEL 
1. Original Design 

Experiments carried out by Sarpkaya and Tuter [100] on small 
smooth cylinders at low Reynolds numbers, in 1974, showed the necessity 
of an apparatus capable of producing larger Reynolds numbers. Various 
designs were considered, ranging in a broad spectrum of complexity and 
versatility. One of the options was to oscillate the fluid by a suitable 
piston mechanism. The achievement of desired amplitudes and frequencies 
needed to obtain high Reynolds and Keulegan-Carpenter numbers required an 
extremely complex system. Furthermore, the vibrations generated by the 
piston motion could not have been easily overcome. Thus, this design has 
been abandoned. Finally, among various options, a U-shaped oscillating- 
flow tunnel was designed and built [101]. The tunnel in its original 
design consisted of 1] modules for ease of construction, transportation, 
and mounting. Each module is made of 3/8 inch aluminum plates and rein- 
forced with 1/2x4x18 inch aluminum flanges welded to the plates. The 
modules were assembled with the help of an air dyring silicon-rubber seal 
between the flanges of two adjacent modules and | inch steel bolts placed 
6 inches apart. The inside of each module was precision machined so that 
the largest misalignment was about 0.04 inches. 

The cross-section of the two legs is 6x3 feet whereas that of the 
test section is 3x3 feet. This selection was dictated by several factors 
Such as: 

(a) available ceiling height; 


(b) pressures to be encountered; 
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(c) desire to obtain, at the test section, an amplitude or 
velocity of oscillation at least twice that of the free 
Surface; 

) period of oscillation; 

) desired Reynolds numbers and relative amplitudes; 

(f) possible separation of the flow at the corners; 

) natural damping of the oscillation; and 

) magnitude and the frequency of the forces. 

The length of the horizontal test section was chosen larger than twice 
the actual amplitude (amplitude at the test section or virtual amplitude) 
to ensure fully developed uniform flow at the test section. Finally, the 
two corners of the tunnel were carefully streamlined to prevent flow 
separation (see Fig. 4). This design proved to be more than adequate for 
the experiments planned at that time. 

The auxiliary components of the tunnel consisted of plumbing for 
the filling and emptying of the tunnel with hot and cold water (50° to 
120°F), a heat exchanger, butterfly-valve system, and the air supply 
system. The tunnel capacity is about 6000 gallons and could be filled 
by city water or by means of a pump from a reservoir. 

The butterfly-valve system (mounted on top of one of the legs of 
the tunnel) consisted of four plates, each 18 inches wide and 36 inches 
long. A 1 inch steel shaft was placed at the axis of each valve plate 
(see Fig. 5). Aluminum housings supported both ends of the shaft with 
self-aligning ball bearings. A 6 inch gear was attached to one end of 
each shaft which extended beyond the bearing. All four valve plates were 


then aligned and driven by a simple rack and piston system. The rack 
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Fig. 5 Photograph of butterfly valve system 





as actuated by an air-driven piston with the help of a three-way valve 
connected to the laboratory air supply system. 

To set the fluid in harmonic oscillation, butterfly valves at 
the top of one of the vertical legs were closed allowing the sealed 
chamber to be charged with high pressure air. This above-ambient pressure 
caused the water tunnel to rise in the opposite leg, left open to the 
atmosphere. The three-way control valve enabled one to release the air 
quickly by sudden opening of the butterfly-valve and put the bulk of 
fluid into a harmonic oscillation. The damping of the fluid oscillation 
was in a very acceptable range in the case for the experiments with the 
fixed cylinders. But, as will be seen later, design improvement was 
necessary in the case of oscillating cylinders. 

2. First Design Improvement 

In order to investigate the effect of cylinder oscillations on 
the transverse forces acting on it, it was necessary to construct an 
elastic mounting with only one degree of freedom. Several design options 
were considered. Finally, the requirements for rigidity, strength and 
torsional stability in the in-line direction quickly led to the use of 
helical steel springs. In the final design, the cylinder was rigidly 
fixed to two low-friction linear bearings resting on two matched helical 
Springs. The motion of the bearings was constrained by two rods perpen- 
dicular to the direction of flow as shown in Fig. 6. 

The force transducers were mounted as shown in Fig. 7 to the 
aluminum split-bearing housing blocks. Each block was manufactured to 
accommodate an adjustable ball bearing with a nominal 2.5 in. bore and a 


length of 5.0 in. The bearings rode on 2.5 in. diameter stainless steel, 
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Fig. 6 Preliminary sketch of the support and housings 
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Fig. 7 Side cutaway view of a load cell and bearing housing, 





hardened and ground, solid shafts. Adjustments of the bearing diameter 
eliminated bearing chatter. The springs were centered around the shafts 
by a retaining ring on the housing blocks and by the yoke support at the 
bottom. 

The yoke supports allowed the blocks and springs to be raised or 
lowered to center the cylinders vertically in the test dation and per- 
mitted the use of various cylinders, springs and weights while accommo- 
dating the effects of the buoyant force. Each yoke was positioned by 
another yoke which was positioned with precision by a threaded third brass 
rod passing through the upper yoke and lifting against the housing top. 
Alignment of mounts and rods was accomplished by a jig which was placed 
between the two blocks, leveled, then raised, and lowered to position the 
two stainless steel rods. Detailed drawings of the system are shown in 
Figs. 8 and 9. 

Figure 10 shows a dimensioned drawing of one of the two support 
housings. Each was welded of 3/8 in. aluminum plate. The windows were 
covered by removable Plexiglas plates for access for calibration and 
maintenance, or, when mounted in place, for the observation of cylinder 
movement during tunnel operation. Vents and drains were installed at 
top and bottom of the housing to remove air or water in filling or draining. 
The housings were flanged and bolted to the U-tunnel wall with O-ring 
Seals for water proofing as may be seen in Fig. 9. 

A 6.0 in. circular access hole was cut in the tunnel walls at 
the test section. Each hole was covered with a plate of identical size. 
A slot of 1 in. by 6 in. was cut in each cover. The covers allowed the 
installation of the cylinders; and the slots, the unobstructed linear 


oscillation of the force-transducer arms attached to the cylinder. 
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Fig. 8 Side view of a cylinder-support mount and housing, 


(dimensions in cm) 
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Fig. 9 Top view of a cylinder-support mount and housing, (dimensions in cm) 
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Fig. 10 Support housing, (dimensions in cm) 
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The usefulness of this design stems from the fact that the 
support-housing assembly could be rotated 90 degrees, thus enabling one to 
conduct experiments for cylinder oscillations in the in-line direction. 

3. Second Design Improvement 

A body oscillating in a fluid will extract energy from it. 
Therefore one might expect this energy loss would add up to the natural 
damping of fluid oscillation in the tunnel giving rise to a faster rate 
of decay of fluid oscillation amplitudes. This would mean that lock-in 
phenomenon could not be observed over many cycles of oscillation. Experi- 
ments performed by Klinkhamer [102] suggested the necessity of devising a 
system capable of restoring the energy lost, due to damping, to the fluid. 
This would enable one to maintain the amplitude of the flow oscillation 
at a desired level thus allowing the observation of the lock-in over a 
large number of flow cycles. For this purpose an air-supply system was 
designed to introduce large quantities of air in one of the legs of the 
tunnel during prescribed time intervals. The necessary air was supplied 
by a centrifugal fan driven by a 3HP electric motor. The fan was first 
connected to a 30x30x33 inch air-control box which regulated the quantity 
of air supplied to the tunnel. The air-control box consisted of a gate 
valve at the exhaust of the fan and another gate valve on top of the box. 
The first gate valve was raised and lowered vertically by a threaded rod 
attached to a manually operated gear. This allowed a fine adjustment of 
the air supply. The gate valve on top of the box was activated by moving 
it horizontally by hand for a coarse adjustment of the air supply. 

The air was introduced into the tunnel by a duct through a simple 
control system (see Fig. 11). The function of the control system was to 


introduce the air into the tunnel at prescribed time intervals so that 
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Fig. 11 Sketch of the air Supply system 
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the energy lost due to damping could be restored to the fluid (see Pig. 
12). The air control system consisted of a drum of 2 feet in diameter 
and 3 feet long. Three openings at angles as specified in Fig. 12 were 
cut into the drum. The door number (1) functioned as inlet to the tunnel 
and door number (2) worked as exhaust to the ambient. Inside the drum 

a vane of the same size as the diameter and length of the drum was placed. 
The vane was rocked about a mean position by an adjustable connecting 

rod through a sprocket wheel and chain mechanism. The system was driven 
by a variable speed electric motor and gear. Two flywheels, one attached 
to the drive shaft of the connecting rod and the other to the motor shaft 
assured smooth and uniform rotation. An RPM counter switch was installed 
at the other end of the drive shaft and the signal was electrically trans- 
mitted to a digital counter on the instrument panel. 

The top of the leg of the tunnel to which air was supplied had to 
be sealed to prevent air from escaping. This was done by plywood boards 
in such a way so as to create an air-tight chamber on the top of the 
tunnel leg. 

With this system very high amplitudes of flow oscillation were 
achieved. This, in turn, necessitated the lengthening of the horizontal 
Section of the tunnel in order to maintain uniform flow about the test 
cylinders. For this purpose two additional modules (each 2.5 ft long) of 
the same cross section as the existing test section were added to each 
Side of the test section. On both sides of these new modules two side 
windows (14x19 in.) were cut. One of the modules also had a door on top 
of it. The side windows were covered with plexiglass for the observation 


of the cylinder motion. However the main purpose of these windows was 
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to provide an easy and quick access into the tunnel for the mounting of 
the cylinders and force calibrations. Finally, a new water level sight 
window was cut into the wall of the vertical leg to facilitate the reading 
_ of the water level. The old butterfly-valve system on top of the other 
leg was not removed since by closing it the flow oscillation could be 
quickly stopped. 

Figures 13a and 13b show photographs of the two sides of the tunnel 
after its completion. Figure 14 shows a schematic drawing of the tunnel. 

The natural period of the tunnel after these modifications turned 
out to be 5.960 seconds. 

To put the fluid in the tunnel into oscillation, the RPM of the 
variable speed motor driving the air-valve was set at such a value that 
the period of the rotation of the vane was exactly the same as that of 
the tunnel. Once the fan was turned on, every time the vane cleared 
port (1) of the air-valve, air flowed into the tunnel and pushed the water 
down. The time of opening of the port (1) was exactly synchronized to 
the time when the water had reached its maximum amplitude and began to 
fall down. The pressure of the supplied air restored the necessary 
energy to the water to maintain its amplitude at a certain level. By 
Opening or closing the gate-valves on the air-control box, one could 
change the quantity of air and consequently the amplitude. This design 
turned out to be extremely successful because the tunnel could sustain 


oscillations at a desired amplitude for an indefinite period of time. 


B. CIRCULAR CYLINDER MODELS 
Circular cylinders with diameters ranging in size from 2.5 in. to 


5 in. were used in the experiments. The cylinders were turned on a lathe 
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Fig. 13a Completed tunnel with air supply system 
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Fig. 13b Back view of the tunnel with the side windows 
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from aluminum pipes and polished to a mirror-shine surface. The length 
of each cylinder was such as to allow a 1/32 in. clearance between the 
tunnel wall and each end. As part of the total length, smooth circular 
flanges 6.0 in. in diameter and 1.8 in. thick were attached to the ends 
of the cylinder to eliminate the end effects. 

Cylinders were hollow, and each contained a beam type accelerometer. 
Figure 15 shows a cutaway drawing of a cylinder with the accelerometer 
installed in the exact center of the cylinder. The electrical connection 
to the accelerometer was led. through a notch in the male fitting at the 
end of the cylinder. Sealing of the unit was accomplished by O-rings 
| and threaded fittings. 

The same cylinders were also used as rough cylinders by coating them 
with sand of desired size. The relative roughness used in the present 
experiments was k. =e . For this purpose clean sand was screened 
through the use of a combination of appropriate size of sieves to separate 
the desired size. The cylinder surface was then cleaned with alcohol and 
carefully coated with an epoxy paint. It was possible to obtain an 
extremely thin and uniform coating which was free of ridges, waves and 
local buildup. Next,. the selected sand was sprinkled over the cylinder. 
The excess sand was easily removed after 15 minutes through gentle rubbing 
by hand. The resulting roughness was uniform as evidenced by the two 
sample photographs shown in Figure 16. The pictures of the rough suface 
were taken with a scanning electron microscope. The actual size of the 


Sample sand is shown in the figure. 


C. FORCE MEASUREMENTS 
Two identical force transducers, one at each end of the cylinder, 


were used to measure the instantaneous force. The basic transducer was 
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Fig. 16 Scanning Electron Microscope Photographs of sand-roughened 
surface 
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manufactured by B.L.H. Electronics, Inc., under the trade name of LBP-1. 
Two sets of transducers were used in the experiments. Each gage in one 
of the sets had a capacity of 500 1b and in the other set a capacity of 
250 1b. The overload capacity of the gages was 200 percent. The deflec- 
| Hon of the gage under 500 Ib load was 0.01 inch. A photograph of the 
basic transducer is shown in Fig. 17. The bellows which protected the 
strain gages had to be waterproofed in such a manner that they would not 
adversely affect the operation of the gage when subjected to about 20 
feet water pressure at temperatures ranging from 65°F to 160°F. For this 
_ purpose the bellows were completely filled with Dow Corning 340-RTV 
coating without bringing the rubber into contact with air during the 
filling operation. After filling, the ends of the bellows were sealed 
air tight with special clamps. The transducers then were mounted, as 
shown in Fig. 7, to the aluminum split-bearing housing blocks. 

It is seen from Fig. 18 that the transverse force acting on the 


cylinder is given by 


in which MX represents the inertial force due to the acceleration of 
the cylinder. The sensitivity of the accelerometer was adjusted to yield 
FE = 0 while the cylinder freely oscillated in air. This method allowed 


the measurement of the net fluid force with an error less than 1.5 percent. 


D. CYLINDER DISPLACEMENT MEASUREMENTS 
To measure the transverse displacement of the cylinder during its 


oscillations a linear-variable-displacement transducer (LVDT) was installed 
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Fig. 18 Free-body diagram of the test cylinder 
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between each bearing block and housing. Each LVDT was energized by a 

six volt DC regulated power supply. The LVDT's provided displacement 

from the equilibrium position and also helped to verify the simultaneous 

and in-phase oscillation of both ends of the cylinder (see Figs. 19a through 


19c for sample displacement, force, and acceleration traces. 


E. ACCELERATION, ELEVATION, OR VELOCITY MEASUREMENTS 

It is because of the extreme importance of the accurate measurement 
of the instantaneous value of these quantities that they are discussed 
here separately. 

It should be noted that the measurement of the amplitude, accelera- 
tion, elevation, or the velocity is a matter of interpretation of the 
Signal received from the appropriate transducer in light of one of the 


following expressions. 


S: >. i? a" ya = Bu, (3) 
in which T is the tunnel period and is a constant. In these experiments 
use was made of a differential pressure transducer with pressure taps 
located symmetrically on the two legs of the tunnel at an elevation 50 
inches below the mean water level, i.e., H = 38 in. Applying Bernoulli's 
equation for unsteady flow between each pressure tap and the instantaneous 
level of water, it is easy to show that twice the amplitude of the free 


surface oscillation (virtual amplitude) is given by 
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Fig. 19 Sample displacement, force, and acceleration traces 





in which g and T are constant and H is kept constant (see Fig. 14). 
Thus the signal of this transducer yielded the virtual amplitude or 
maximum velocity in each cycle. It was entirely free from noise or 
free surface effects. The transducer was calibrated and its linearity 


checked before each series of experiments. 


F. MEASUREMENT OF OTHER BASIC PARAMETERS 
The total stiffness of the system, Kos was determined by suspending 
a known weight from the center of the cylinder, and the displacement 


was measured. Then Se was calculated from the relation 


K Force (5) 


s Displacement 


The natural frequency in air, f.., and the natural frequency in 


na 
water, ae were determined by plucking the cylinder and allowing it to 
oscillate freely in each particular fluid medium. The displacement was 
recorded on one channel of a three channel recorder. A sample trace is 
shown in Fig. 19d. Along with the frequencies, the logarithmic decrement 
and the damping factor, Ca? and G42 were determined in air and water, 
respectively. 

Total mass of the system was measured by adding up the mass of the 


components of the oscillating system. The system mass was also calculated 


from 
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Fig. 19d Displacement trace of a "plucked" cylinder. 





to check the correctness of the frequency evaluation. The results were 
identical. 


Added mass We of the system was determined from 


S 
M, +M = = M es) (7) 
(2nf_ ) 
nw 
For the purpose of verification of the measured values, Mo for the 
cylinder was computed through the use of the observed values of frequencies 
and found to be within one tenth of one percent of that given by the 


potential flow theory and independent experiments. 


G. INSTRUMENTATION AND ELECTRONIC CIRCUITRY 

The output of each of the two LVDT's, after amplification by a 
low-gain amplifier, was displayed on an eight-channel recorder as in 
Fig. 20. One of the two signals was additionally displayed on a three- 
channel recorder. | 

The output of the force transducers was separately amplified by a 
carrier amplifier, filtered by a Krohn-Hite electronic filter and dis- 
played on an eight-channel recorder. Then the two signals were added 
together by an electronic summing circuit, amplified again by a low 
level preamplifier and displayed on the eight-channel recorder. 

The accelerometer output was similarly amplified by a carrier 
amplifier and filtered by a low pass filter. This signal was also 
displayed on the eight-channel recorder. It was then subtracted from 
the sum of the force signals through the use of a summing-differencing 
electronic circuit. The resultant signal, representing the fluid force 


on the cylinder, Fis was again amplified by a low level preamplifier and 
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displayed on both three- and eight-channel recorders. Figure 21 is a 
block diagram of the circuit described above. 

The signal from the differential pressure transducer was of such 
quality that it did not require any filtering. As shown in Fig. 22, 
jit was amplified by a carrier amplifier and displayed on the two afore- 


mentioned recorders. 


H. PROCEDURES 

Five cylinders were used in these experiments. Additionally two 
Springs, manufactured to the specification desired, were used. The 
diameter of the springs were sufficiently large to work freely over the 
2.5 in. diameter rod. The coil diameter of the springs was 0.255 inches 
and were made of carbon steel wire. Table 3 is a condensed list of the 
cylinders used, kind of force measured, surface condition, and spring 
constant employed in the experiment. The number of (x) in each box 
indicates the number of experiments performed under that particular 
conditions. The check mark under the spring-constant box indicates the 


kind of spring used. 


Table 3 
Diam. LIFT DRAG K. (Lb/ft) 
in. Smooth Rough Smooth Rough 418 392 
a X v 
2.5 K X vy 
B. XXX i X X vy 
4 ro. KK / / 
5 OK XX / 
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Fig. 22 Differential pressure transducer and circuitry 
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In some of the experiments it was found advantageous to increase 
the weight of the oscillating system by mounting rectangular lead blocks 
to the side of the oscillating bearing blocks away from the cylinder 
mountings. Beside adding mass and therefore lowering the natural fre- 
quency, it was discovered that adding the weight at this location also 
served to distribute the load better on the bearings, thereby reducing 
friction and damping. 

Procedures for an experiment were as follows: 

1. The cylinder was aligned with the tunnel by adjusting the 
bearings. The sensitivities of the force transducers were matched by 
hanging a known weight from the exact center of the cylinder. At the 
same time a calibration factor was determined by measuring the displacement 
and comparing it with the output of the LVDT. 

2. The accelerometer was calibrated by plucking the cylinder and 
adjusting the amplifier gain until the fluid force in air was nulled to 
within 1.5 percent deviation. 

S. fey and Gc, were also determined by plucking the cylinder in air 
from the trace of cylinder oscillation on the three-channel recorder 
(see Fig, 19d). 

4. The weight was then suspended from the cylinder center, this 
time for the measurement of a calibration factor in 1b/mm of the three- 
channel-recorder output at a particular gain setting. 

5. The U-tunnel and cylinder-support housing were again filled to 
the top of the test section, and all air was vented from the housings. 

6. The cylinder was then lowered to the center of the test section 


by simultaneously lowering both sets of spring yokes, each the same 
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number of turns, to compensate for the buoyant force on the displaced 
cylinder. 

7. The cylinder was plucked, and ee and c, were calculated from 
the free-oscillation traces. 

8. The access to the horizontal test section above the cylinder, 
through which the weights were lowered to the cylinder, was closed. 
The tunnel was filled with water from a sump tank located near the 
apparatus or from the city water. 


9. The differential-pressure transducer was bled of any entrained 


’ air and its calibration verified. 


10. Variable speed motor control, driving the air valve, was adjusted 
to yield a period exactly equal to that of the tunnel. 

11. Butteryfly-valve on top of the air chamber and the gate-valve 
On top of the air control box were fully closed. The gate valve at the 
outlet of the fan was fully opened, and the butterfly-valve on top of 
the other leg was fully opened before the fan was turned on. With this 
arrangement the largest amplitude of the tunnel oscillation was obtained. 

12. Several tunnel oscillations were allowed to pass in order to 
achieve a steady state tunnel oscillation. 

13. Amplifier gain settings or calibration factors were recorded 
On the chart paper as they were changed as desired. 

14. Several cycles of oscillation were recorded at a chart speed 
of 2 and 10 mm/sec to ensure that a steady state is reached for that 
particular amplitude of the tunnel. 

15. Several runs were made at a chart speed of 25 mm/sec. 

16. The position of the gate valves on the air control box was 


changed to achieve a new amplitude of the tunnel oscillation. 
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17. Steps 14 through 16 were repeated. 


18. At the end of each experiment, Lo oy? f _, and G4 were measured 


na 
to verify the initial values. 

In measuring the forces in the in-line direction, steps 1] through 
4 were repeated. The force transducers were turned 90 degrees and pinned 


to measure the in-line force. Then the steps 5 through 17 were followed 


as before. 
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22.2 








ITT. METHOD OF DATA ANALYSIS 


A. INTRODUCTION 

The recorded data consisted of the instantaneous values of the 
transverse displacement of the cylinder, lift force, and the acceleration 
of flow. Figures 19a through 19c¢ show a sample of the traces of the above 
variables during twe cycles of flow oscillation. 

The data were subjected to extensive analysis in terms of the peak 
amplitudes of the force and cylinder displacement, their root-mean-square 
values, Fourier analysis, spectral analysis, Morison's equation [96], 


and finally, in terms of Duhamel's superposition integral. 


B. ANALYSIS OF THE PEAK AMPLITUDES 


The parameters D, Ks aa re Ca? and Cy were recorded for each 


run. Then the peak values of the transverse force (lift force), Fa? 
peak values of the cylinder displacement, Ku and the periods Te and 
T. for the force and displacement, respectively, were evaluated over each 


half cycle of flow oscillation. Subsequently, the variables A. UT 


K, Xy/ Ds ULF Ds and C,,, were calculated for further study. The coefficient 


LM 
CM is the peak lift coefficient and is given by 


; LM | 
cM TD (8) 


7eU Louk 


where o is the density of the fluid; and L, the length of the cylinder. 
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C. DIGITIZATION OF THE DATA 

The evaluation of the root-mean-square values of the various 
parameters, Fourier analysis, spectral analysis, Duhamel's superposition 
integral, and the use of Morison's equation required the careful digiti- 
zation of the data. 

In order to maintain a proper relationship between the phase and 
amplitude of a given variable a suitable frequency of digitization or 
digitization interval had to be chosen. Otherwise aliasing could occur 
which is a potential source of error in an analog to digital data conver- 
sion [103]. The maximum frequency resolvable is called the Nyquist 


frequency and.is defined as 
7 oe (9) 
NY 2At 


where At is the interval of digitization. The frequencies in the original 

data above this cut off frequency are folded back into the frequency range 

from zero to Fay in an accordion pleated fashion. Therefore, in order 

to choose a suitable At the maximum frequency encountered in the experi- 

ments had to be known. The observation of the analog records showed 

that the highest frequencies encountered were in the range of 4 to 6 Hz. 
In general, it is a good rule to select Fay to be one-and-a-half or 

two times greater than the maximum anticipated frequency. Therefore, a 

Nyquist frequency of 10 Hz was chosen from which an interval of digitiza- 

tion of 0.05 seconds resulted. This corresponded to a 1.25 mm interval 

On the chart at its speed of 25 mm/sec. Any digitization interval below 

this value would have given correct results. The smallness of this 


interval is defined by the economy of the effort. 
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In case of the use of data in a spectral analysis the number of 


digitized data had to be N which is defined as 


N = 2x 2? ee a ee 


This consideration and others such as the ease of tapering led to an 
interval of digitization of 1.073 mm, equivalent to At = 0.043 seconds. 
This time interval is lower than the conservative value of 0.05 seconds. 
The length of the records were different but never less than three 
cycles of tunnel oscillation which is about 18 seconds. This resulted 
in a resolution of 1/18 = 0.056 Hz. The digitized records were long 
enough to satisfy all the assumptions inherent in the analysis. 
The analog traces at a fixed tunnel amplitude were digitized by 
means of a Hewlett-Packard digitizer model 9874A which was linked to a 
Hewlett-Packard 9845A computer. The movement of a cursor by hand over 
the trace generated digitized coordinates of a given curve. The informa- 
tion was then stored either in files on tapes for use with the programs 
written for the HP computer or punched on cards to be used with an I.B.M. 


360 computer. 


D. EVALUATION OF THE ROOT-MEAN-SQUARE VALUES 
The use of digitized data together with a suitable computer program 
yielded the rms values of the lift eoefficient, (Cr mg? and the relative 


cylinder displacement, (X) nef D- 
The record length for the evaluation of the rms values was at least 
three cycles of tunnel oscillation. A careful examination of the data 


has shown that a record length of three cycles was more than sufficient. 
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E. SPECTRAL ANALYSIS 

As noted earlier, the nature of the lift force and cylinder-displace- 
_ment traces were quite regular, particularly in the synchronization range. 
The quasi-harmonic shape of the said traces required a deterministic 
rather than random analysis. However, a knowledge of the frequency ccn- 
tent of the oscillations was of primary importance in determining the 
modulations in the amplitude and frequency of the oscillation of the 
lift force and cylinder displacement. Evidently, only through spectral 
analysis can one decompose the force and displacement into their funda- 
mental and subharmonics and draw conclusions regarding the energy content 
of each harmonic. 

The following is a condensed description of the methodology. 

If x(t) is a periodic function of time t, with a period Ts then 


ijt could always be expressed as an infinite trigonometric series of the 


form 
- 2nnt __,2nnt 
x(t) = a, + ) [a_ cos( = | eee sin) (10) 
n= S S 
where 
- 
Lf” x(t) dt (11) 
a= = x 
0 lie 5 
I 
a = a alc), cos (27) dt (12) 
a3] SO S 
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b= - f x(t) sin(Stmty at (13) 


n>] 


The mathematical conditions for the convergence of Eq. (10) are extremely 
general and cover practically every conceivable engineering situation 
[104]. 

The Fourier series is a useful tool for determining the frequency 
content of a time-varying Signal. However, the Fourier series always 
require a periodic time function. To overcome this shortcoming, Fourier 
evaluated his series as he let the period of the waveform approach 
infinity. Then the Fourier series given by Eq. (10) turns into a Fourier 
integral and Fourier coefficients represented by Eqs. (11), (12) and 
(13) become continuous functions of frequency, called Fourier transforms. 
Therefore, a Fourier integral may be regarded as the formal limit of a 
Fourier series as the period tends to infinity. It is customary to 
write the Fourier series in their complex form. Then one introduces 


the Fourier transforms defined by 


FL (f) = met) ge léntt dt (Forward Fourier transform) (14) 
x(t) = f at) eléntt df (Inverse Fourier transform) (15) 
where 
enucunt = cos(2nft) + i sin(2nft) , 
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is known as the Kernel of the Fourier transform. Now FCF) contains 

the amplitude and phase information at every frequency present in x(t) 
without demanding that x(t) be periodic. If the time history x(t) of 
a sample function is not periodic it cannot be represented by a discrete 
Fourier series. Furthermore, for a stationary process, x(t) goes on 
forever and the necessary condition 


co 


f |x(t)| dt < @ (16) 
is not satisfied. Unless special precautions are taken, it is not 
therefore possible to evaluate the Fourier transform of x(t) in order 
to obtain information about the frequency composition of the process. 
This difficulty can be overcome by Fourier analysing, not the sample 
function of the process itself, but rather its autocorrelation function 
R(t). 

The logic behind this approach is that the autocorrelation function 
gives information about the frequencies present in a random process 
indirectly. If the zero value of the process x(t) is normalized (or 
adjusted) so that the mean value of the process is zero and it has no 


periodicity, then the autocorrelation function satisfies Eq. (16) and 


its Fourier transform could be written as 


ee) 


mes ff Rit) ee dr (17) 


=O 


The S (fF) is called mean square spectral density, energy density, or 


power spectrum. The last name originated from electrical engineering 
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because it is a measure of the power loss ina resistor. It can be 


shown that S.(f) is a real, even function of frequency, f and it is 


always positive [103, 105]. 


One sided spectral density is defined as 


G(f) = 2S (f) 


It can be shown that 


where a, @ 1s the variance of the process [106]. 


The one sided cross spectrum is found from 


(f) 


G 
xy et xy 


or 


mee) = ae (Tt) 


xy xy a, Qey 


where Cy 


The Eat) is the quad-spectrum and forms the imaginary part or out of 


phase component of G yi Ff). 


From Eq. (21), the phase difference between the process x and y is 


found to be 


ae 


2 f R. Ge IRA 


(f) 


(21) 


(f) is the co-spectrum. It is the real and in phase component. 





(fF) 


ia 
' AL. (22) 


Oy = tan C 


X 
XY 


If a finite record of length ie is used, it can be shown [103] that an 


‘estimate of the cross-spectral density is given by 


Best) = 2 nL in 


Xy S AGP (23) 


where X and Y represent the finite Fourier transforms of x and y respectively 
and X* is the complex conjugate of X. Then an estimate for spectral 


density of a Stationary process of finite length would be 


2 
G(f.T.) = 2 an ; (24) 
where 
Xin) = x(a, - ib ) (25) 


in which ay and b are the Fourier coefficients given by Eqs. (11), 

(12), and (13). If N digitized data are used to evalute ay and b, 
numerically then the computation time is proportional to Ne. Even with 
high speed computers, computation of discrete Fourier transforms requires 
excessive machine time for large N. In 1965 Cooley and Tukey published 
their mathematical algorithm [107] which has become known as the "fast 
Fourier transform". The fast Fourier transform (F.F.T.) is a computational 
algorithm which reduces the computing time of ay, and b, to a time 


proportional to N log N. The idea of F.F.T. was first discovered by 
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Runge [108, 109] in 1903 who described essentially the computation of 
the F.F.T. algorithm as we know it today. 

Any spectral estimate has only two degrees of freedom and the 
confidence interval turns out to be so large that it does not have any 
statistical meaning. Therefore, one deduces that the evaluation of 
only one G (fT) for a record length does not have any significance. 
There are several methods to increase the number of degrees of freedom. 
A simple method is to create an ensemble of records each of which has a 
record length T.. If the length of the original record is long enough 
one can divide the record into k records each be seconds long. Now 
ensemble averaging improves the degrees of freedom, which in turn improves 


confidence, or spectral reliability of the spectral estimate then 
rs — 
G(f.T.) =~ } G(f,T..k) (26) 


where k is the number of sample function forming the ensemble and the 
degrees of freedom for each spectral estimate is now A = 2k. The true 


spectrum is found as 


kK A 
) G(f,T_sk) (27) 


G(f) = Lim 
T 


s k 
k->00 


The use of this method requires that the process be ergodic. This means 
that the process must be stationary, i.e., for any sample function k, 
time average of the mean, autocorrelation and all other statistical 


averages do not differ when computed over different sample functions. 


94 





As it was mentioned earlier, the observation of the force and cylinder 
displacement traces revealed that the phenomenon repeated itself at 
almost every three cycles of flow oscillation for a given Keulegan- 
Carpenter number. This justified the assumption of ergodicity for the 
process. 

Two or more sample functions each of which comprising at least 
three cycles of tunnel oscillations were used to evaluate spectral 
density of the process at a constant Keulegan-Carpenter number. 

A flow chart of the program used to evaluate the spectral density 
is given in Fig. 23. The data was first "detrended" to remove any linear 
trend from the process. Then a cosine taper window was used to reduce the 
leakage and smoothen discontinuities at the "joints" at the two ends of 
the records. 

In addition to the spectral density for each data set, other statisti- 
cal averages such as variance, standard deviation, skewness, and kurtosis 
were estimated. Integration of spectral density, using Eq. (19), was 
performed to compare the result against variance which was evaluated 
directly from digitized data. The closeness of these two was also a 


good check for the goodness of spectral estimate. 


m IN LINE FORCE 

The drag coefficient Cy and the inertia coefficient C. for the in- 
line force were evaluated in order to determine the effect of the cross- 
flow oscillation of the cylinder on the in-line force. Data reduction 
for the forces in-line with the direction of the ambient flow was based 
on Morison's equation [96] and the Fourier analysis. The in-line force, 


consisting of the drag force F and the inertia force Fs is assumed 
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Fig. 23. Flow chart of the computer program for the 
evaluation of the spectral density. 


96 











to be given by 


; i i 2 du 
mete ao, L Do Ul) U+ee LO nos (28) 


in which Cy and Ces represent respectively the drag and inertia coeffi- 
cients and U is the instantaneous velocity of the ambient flow. For an 


harmonic flow represented by 
U = -U, cos @ (29) 


where 


equation (28) may be written as 


2 
Tl : 
= -C, |cos @| cos 8+ C ~ sin ® (30) 


The Fourier averages of Cy and Cf were evaluated from (30) by Keulegan 


and Carpenter as [I11] 


3 oT Re 
Cy ae f ——5z—— cos 8 dé (31) 
0) 0 U DL 
° m 
So ail F 
c = Sf N___ sin @ do ea) 
i O fe) U bee 


in which be represents the measured force. 
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The in-line and lift forces were measured by two force transducers 
at each end of the cylinder while it was oscillating in the transverse 
direction. The three channel recorder with three strip chart paper was 
used to record forces and the flow oscillation. Figures 24 and 25 show 
sample in-line and transverse force traces. The in-line force was then 
digitized every 0.1 seconds and the digitized data, along with other 
parameters such as D, L, vy, and A were fed into a computer program to 


evaluate Cy and C. from Eqs. (31) and (32). 


G. DUHAMEL 'S SUPERPOSITION INTEGRAL 

One of the methods used to investigate the relationship between the 
instantaneous force and displacement in a linear, oscillating, system 
is the so-called Duhamel's superposition integral. Assume that the 


Oscillation of the cylinder is represented by 
MX+2Mc 0 K+k X = F 


where F(t) is the instantaneous lift force. One can imagine that the 
force F(t) is composed of a series of impulses [111]. Then through 
the use of the convolution or superposition integral one can evaluate 


the instantaneous response of the system to the force F(t) from [111] 


lee, b(t) ear fac 
_ eC ] - 
= (X, rae f sinyl-z wt dt) cos wt 
n oO ime z 

-tw (t-t, ) 

X ic n V 

o, | F(t) e We 

+ a re f : cos ¥1-z wt dt, )sinw t (34) 
n no | oe 
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Fig. 24 A sample of the in-line force trace. 
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Fig. 25 A sample of the transverse force trace. 
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in which Xo and - denote the initial condition and CY is the variabie 
time. 

For this purpose, the lift force trace was digitized every .02 
seconds and used as the input to Eq. (34) along with other parameters 
characterizing the oscillating system such as cz, D, M, etc. A computer 
program then performed the calculations. This computation served also 


as a confirmation for the consistency of the lift force data. 
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IV RESULTS 


A. INTRODUCTION 

The experimental results will now be reported, showing first the 
dependence of the relative displacement and lift coefficient on the 
flow parameters, and, secondly, the analysis of the data through the use 
of the Fourier analysis, spectral analysis, and Duhamel's integral. 
Lastly, the in-line force acting on transversely oscillating cylinders 


will be presented. 


B. REPRESENTATIVE RELATIVE DISPLACEMENTS AND LIFT COEFFICIENTS 

The variation of the relative displacement Xy/D as a function of 
K = U, iD) U. = ULF Ds Lie le and Val Nie 1s shown in Figs. 26 through 
29 for a 3 in. smooth cylinder. Data for other representative smooth 
and rough cylinders are shown in Figs. 30 through 45 in Appendix A. In 
the foregoing, U., represents the maximum velocity of flow during a cycle; 
Fes the frequency of the transverse-force oscillations between two 
Successive peaks; os the natural frequency of the cylinder in still 
water; and Tae the frequency of the cylinder oscillations. The relative 
roughness for all rough cylinders was kK /D = 1/100 where Kk. is the 
average sand-roughness height. 

The variation of the rms values of X/D with K, Us TAL sis and Tho. 
for the 3 in. smooth cylinder is shown in Figs. 46 through 49. Similar 
data for other cylinders are shown in Figs. 50 through 65 in Appendix B. 
The peak amplitudes of the lift coefficient are shown in Figs. 66 through 


69 as a function of K, Us Ved aaa and F/T. for the 3 in. smooth cylinder. 
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Data for other cylinders are shown in Figs. 70 through 85 in Appendix C. 

The rms values of C. for the 3 in. cylinder are shown in Figs. 86 through 
89 for the 3 in. cylinder. The remainder is shown in Figs. 90 through 
105 in Appendix D. 

The data presented in the figures cited above exhibit certain 
characteristics which will be discussed below. Firstly, all the plots 
Showing the peak-amplitude values of the relative amplitude and the lift 
coefficient show considerable scatter. This is partly due to the amplitude 
and phase modulations in the force and displacement traces and partly due 
to the inherently random nature of the vortex shedding in steady and 
harmonic flow about stationary and oscillating cylinders. The rms values 
of the displacement and force exhibit much smoother variation with respect 
to the parameters K, Un Are and f,/f The data span a broad range 
of Keulegan-Carpenter numbers. Within this range the relative amplitude 
shows one dominant and, at times, one secondary peak. For example, Fig. 
26 shows that Xi/D reaches a maximum value of 0.27 at a Keulegan-Carpenter 
number of about 55. A second and considerably smaller peak occurs at 
about K = 135 (see Figs. 26 and 46). The region in which the first peak 
occurs is called the primary region of synchronization and the region in 
which the secondary peak occurs is called secondary synchronization. 

A tertiary synchronization region observed in steady flow past oscillating 
cylinders has not been found in the harmonically oscillating flow. Clearly 
the primary synchronization region extends over a small range of K values. 
Outside the synchronization region the cylinder continues to oscillate 

at small amplitudes as a body subjected to forced oscillations by shedding 


Of the vortices. The dependence of the primary synchronization region on 
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the reduced velocity Uf td 1S seen clearly in figures similar to 
Figs. 2/ and 4/7. The primary synchronization (referred to hereafter as 
the synchronization) occurs at U. values characteristic of those 
encountered in steady flows [see Fig. 3 where V/ (anf D) = 1]. The 
frequency ratios of f,/f, and f/f (see Figs. 28 and 29 and Figs. 48 
and 49) show that synchronization takes place when the frequency of force 
oscillations nearly coincides with the natural frequency of the cylinder. 
Furthermore, the cylinder oscillates at a frequency close to its natural 
frequency, 

The lift coefficient at synchronization is several times larger 
than that in the non-synchronization region (see, for example, Figs. 
66-69 and Figs. 86-89). As will be noted later, the said lift coefficient 
is 2.75 times that for the stationary cylinder at the corresponding 
Keulegan-Carpenter number. The amplification of the lift force is a 
consequence of the non-linear interaction between the flow and the cylinder 
as discussed in Chapter I. 

As noted earlier, the random nature of the vortex shedding gives 
rise to modulations in both displacement and force traces. It is, 
therefore, useful to examine the relative influence of the harmonics of 


the displacement and lift force. 


C. FOURIER ANALYSIS OF THE SYNCHRONIZATION REGION 

When a signal is periodic over We it is cOnvenient to use a Fourier 
Series representation. The dependence of the Fourier coefficients on 
the pertinent dimensionless groups may then be used to describe the 
Signal completely. Even when the signal is not completely regular, it 


is still useful to investicate the contributions to it, associated with 
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various harmonics of the wave frequency. Thus, if the displacement and 
force signals are considered as given approximately by 


G2) ) xX cos(2nmnft/T + >, ) (35) 


n=] 


and 


ic 


Ln cos(2mnft/T) (36) 


lcd 8 


c(t) = 


n=] 


in which x, and Cy are the effective Fourier components and may be 
deduced through the use of a suitable computer subroutine. 

Figures 106 through 113 show the normalized displacement and their 
corresponding normalized Fourier coefficients for all smooth and rough 
cylinders. An examination of the Fourier-coefficient plot shows that the 
predominant frequency invariably occurs at f = of = 9/T where T is the 
period of flow oscillation in the tunnel. The experiments by Sawaragi 
[99] have also shown that the predominant frequency is equal to a multiple 
of the flow oscillation frequency. The relationship between the predominant 
frequency and Strouhal number will be discussed shortly. Figures 106 
through 113 also show that there are secondary harmonics which are 
considerably smaller than the predominant harmonic. A correct represen- 
tation of the phase and amplitude modulations of the displacement requires 
the use of at least two or three harmonics. Figures 114 through 121 show 
representative normalized force traces and the corresponding normal ized 
Fourier coefficients. 

Once again it is found that the predominant frequency of the lift 


force is f. = 9 Fe Also, there are relatively strong higher harmonics 
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Fig. 106 Normalized displacement and Fourier coefficients 
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107 Normalized displacement and Fourier coefficients 


for 3 in. rough cylinder, ¢ = 
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Fig. 108 Normalized displacement and Fourier coefficients 


for 4 in. smooth cylinder, zc = 0.046. 


124 





ATA POUNTS= (27 
t VHLUE © 1, 996056-01 
Sar rae a | 


a. 0c rt00 


Tatoo) 
ay 
na 


Ag 





NJ TUOE 
ee C§erricrenteo TO 64 


OC TEPNS §. 223536153 
MAX VALUES 1.5%3995E-a1 





| ee = J ae See 


8.9 SS ee ee ere eee 
a 64 
COEFFICIENT 


Fig. 109 Normalized displacement and Fourier coefficients 
for 4 in. rough cylinder, cz = 0.030. 
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Fig.. 110 Normalized displacement and Fourier coefficients 
for 4 in. rough cylinder, ¢ = 0.060. 
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for 5 in. smooth cylinder, zc = 0.067. - 
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Fig. 112 Normalized displacement and Fourier coefficients 
. for 5 in. smooth cylinder, > = 0,055. 
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mages 113 Normalized displacement and Fourier coefficients 
for 5 in. smooth cylinder, ¢ = 0.054, (with lead blocks). 
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Fig. 114 Normalized lift force and Fourier coefficients 
for 3 in. smooth cylinder, z = 0.052. 


130 





1.8 


: ly 


ORTR POLUTS= Lea» 
ABSCHHe VELUC = FS Somers 


-_-—-4 


§.47£168 


e 
ae 





—~ 











=lon .! 


xX VALUE 


1.80 fuss ITI 
VEFFILTENT=H0 75 m4 


oC TEPM= 1.72550 -21 
MAX VALUE= 7. 2oSSE+00 


I ee aoe on oe cere 


COEFFICIENT 


Fig. 115 Normalized lift force and Fourier coefficients 
for 3 in. rough cylinder, ¢ = 0.067. 
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Fig. 116 Normalized lift force and Fourier coefficients 
for 4 in. smooth cylinder, zc = 0.046. 
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Fig. 117 Normalized lift force and Fourier coefficients 


for 4 in. rough cylinder, ¢ = 0.060. 
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5 in. smooth cylinder, z = 0.061, (with lead blocks). 
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on either side of the predominant frequency. It is found that the 
instantaneous lift force, as well as displacement, may be represented 


fairly well through the use of their predominant and two or three 


additional, most significant, harmonics. Figures 122a through 122c show 


the Fourier representation of a sample displacement trace for a 4 in. 
rough cylinder with one, two, and three Fourier components, respectively. 
Similar results for the force trace are shown in Figs. 123a - 123e. 
Evidently, the displacement is fairly well represented by the first 
three components. The force trace for the particular cylinder required 
the use of four Fourier components. 

The occurrence of the predominant frequency at oF is not without 
a reason. If the Strouhal number is defined in terms of the maximum 


velocity Ue then one has 


a (37) 


which varies from about 0.1 to 0.23 during a cycle in harmonic flow. 
Multiplying the numerator and the denominator by T, the period of the 


flow oscillation, one has 








It is fully realized that the vortex shedding frequency Vee does not 
remain constant throughout the cycle. In fact, a is largest near the 
instant of U., and decreases as U approaches zero. Consequently, the 


fundamental question as to which representative vortex shedding frequency 
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the cylinder oscillation synchronizes may be answered only through the 
Fourier analysis. Replacing te by N/T one has SK = N. In the present 
investigation the synchronization occurs in the range of Kelulegan-Carpenter 
numbers from about 48 to 60. With N = 9, one has S = 0.15 - 0.18. 

' Using a value.of S = 0.16 for the perfect synchronization (at which the 
displacement is maximum) one finds K = 56. In other words, the cylinder 
undergoes synchronized oscillations when the natural frequency of the 
cylinder coincides with the vortex shedding frequency based on a Strouhal 
number of 0.16, as verified by the experiments. Consequently the integer 
N may be smaller or larger depending on the Keulegan-Carpenter number 

at which synchronization is realized. For example, for a cylinder whose 
natural frequency is twice the flow oscillation frequency (N = 2), the 
synchronization takes place at K = 12.5. This result is in conformity 
with the results of Sawaragi [99]. However, the Keulegan-Carpenter number 
in his experiments is based on the rms value of the wave induced veloci- 
ties along the vertical cylinder. Thus, an exact verification of the 
corresponding results is not possible. Finally, it should be noted 

that the use of the rms value of the flow velocity and the frequency at 


which perfect synchronization takes place results in S = 0.22. 


D. SPECTRAL ANALYSIS 

The spectral density, coherence, and phase at perfect synchronization 
are shown in Figs. 124 and 125 for the 3 in. smooth and 4 in. rough 
cylinders. Similar results in the non-synchronization region are shown 
in Figs. 126 and 127 for a representative value of K = 80. The spectral 
density for a broad range of K values for the 3 in. cylinder is shown in 
Figs. 128 through 143 in Appendix E. Similar data for the 4 in. rough 


cylinder are shown in Figs. 144 through 158 in Appendix F. 
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Fig. 124 Spectral analysis for a 3 in. smooth cylinder 
at resonance, ~ = 0.052, K = 51. 
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Fig. 125 Spectral analysis for a 4 in. rough cylinder 
at resonance, zt = 0.030, K = 49. 
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Fig. 126 Spectral analysis for a 3 in. smooth cylinder, 
c = 0.052, K = 80. 
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Fig. 127 Spectral analysis for a 4 in. rough cylinder, 
c = 0.030, K = 80. 
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It is clear from Figs. 124 and 125 that the Spectral density is 

- sharply defined at synchronization and the force and displacement peaks 
occur at the same frequency. Furthermore, the coherence is nearly per- 
fect. In the non-synchronized region the spectral density for the force 
becomes broad band and the peaks of the displacement and force do not 
Occur at the same frequency as seen from Figs. 126 and 127. The peak 
of the displacement is relatively better defined and continues to 

occur near the natural frequency of the cylinder fs = 1.69 for the 3 
in. smooth cylinder and fe = 1.54 for the 4 in. rough cylinder). The 
results of the spectral analysis confirm the findings of the Fourier 
analysis that at resonance the force and displacenent maximums occur at 


the same frequency. 


E. DUHAMEL'S INTEGRAL 

Figures 159 through 163 show a comparison of the measured and calcu- 
lated values of the instantaneous displacement as a function of time. 
The measured force shown in the upper half of each figure was used as 
input in Duhamel's integral, as given by Eq. (34). The close agreement 
between the measured and calculated displacements is an indication of the 
proper measurement of the various quantities such as mass, damping ratio, 
frequency, etc. 

A sensitivity analysis was carried out by varying the frequency i 
in Eq. (34) while maintaining all other input parameters unaltered. The 
purpose of this analysis was partly to explore whether the natural fre- 
quency of the oscillating cylinder remained equal to that determined in 
still water and whether the prediction of Duhamel's integral strongly 


depended on the natural frequency. The implication of the equivalence 
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of the natural frequency of the cylinder undergoing transverse oscilla- 
tions in harmonic flow to that in still fluid is that the added mass of 
the synchronized cylinder is identical to the displaced mass of the 
cylinder. If, on the other hand, the natural frequency of the synchronized 
| cylinder were to be higher than that in still water, one would conclude 
that the added mass of the synchronized cylinder is smaller than its 
displaced mass. Extensive calculations with various values of the 
natural frequency have shown that the natural frequency of the synchronized 
cylinder is 5 percent larger than that for the natural frequency in 
still water. The difference between the two frequencies may be partly 
due to the experimental errors and partly due to the difference between 
the actual added mass and displaced mass. King [112] conducted a series 
of experiments in still water (by plucking a cantilevered cylinder) and 
also in flowing water at the fundamental mode of the flow-excited vibra- 
tions. From a comparison of the measured frequencies he concluded that 
the still-water and flow-excited frequencies are virtually equal, ji.e., 
the added mass is unaffected by streaming flow. Apparently, this conclu- 
Sion is valid only at the resonance cOndition. The results of the present 
investigation also showed that the added mass is relatively unaffected 
by the harmonically oscillating flow. One must, however, note that the 
best agreement between the measured and calculated displacement is 
obtained with a natural frequency approximately 5 percent larger than 
the still-water natural frequency. 

The foregoing also points out the important fact that a simple com- 
parison of the natural frequencies measured in still water and in flowing 


water at resonance is not sufficient for an accurate assessment of the 
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Duhamel's integral for a 3 in. smooth cylinder, zc = 0.052. 
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Duhamel 's integral for a 4 in. smooth cylinder, zc = 0.046. 
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Fig. 161 Comparison of measured and calculated displacement by 
Duhamel's integral for a 4in. rough cylinder, ¢ = 0.030. 
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variation or lack of variation of the added mass. The use of Duhamel's 


integral is a more reliable method for the determination of the natural 


7 oe 


frequency at resOnance. Be that as it may, it is rather Surprising that 
the added mass of a body oscillating, at its natural frequency, ina 
fluid otherwise at rest, is nearly the same with that undergoing syn- 


chronized oscillations in a fluid in motion. 


F. DRAG COEFFICIENT FOR TRANSVERSELY OSCILLATING CYLINDERS 
Experiments with cylinders undergoing transverse oscillations in 


steady flow have shown that [113] the in-line force increases with X/D 


_ since the cylinder presents a larger apparent projected area to the mean 


flow. 

From an engineering viewpoint, the significance of the magnitude of 
the in-line force is that a cylinder or cable excited by the flow to 
oscillate in the transverse direction may be subjected to in-line forces 
larger than those assumed in its design. Furthermore, the deflections 
caused by the in-line force of sufficiently flexible cylinders tend to 
couple with transverse oscillations and affect not only the magnitude of 
the transverse oscillations but also the path of the cylinder motion. 

It is because of these reasons that the in-line force acting on the 
cylinders undergoing oscillations in the transverse direction to the 
narmonic flow has been measured. The drag coefficient has been determined 
through the use of Eq. (31) and is shown in Figs. 164 through 166 for 

three representative cylinders. Also shown in these figures are the mean 
lines for the drag coefficient of the corresponding stationary cylinders. 
In general, the oscillating-cylinder drag coefficients are larger than 
those for the stationary cylinder. However the increase is not as dramatic 


as in the case of steady flow about vibrating cylinders [114]. 
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The inertia coefficient C has also been evaluated and it was found 
that it did not differ from that for stationary cylinders. This is 
primarily because of the fact that at Keulegan-Carpenter numbers larger 
than about 30 the in-line force is primarily drag dominated and the 


effect of the inertial force is negligible [see Eq. (30)]. 
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Fig. 164 Cy versus K for a 3 in. smooth cylinder, ~ = 0.052. 
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Fig. 166 Cy versus K for a 5 in. smooth cylinder, x = 0.055. 
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V. ANALYSIS AND DISCUSSION OF RESULTS 


A. INTRODUCTION 

In the following, the results obtained with various analyses and 
models will be compared with those obtained experimentally. First the 
general equation of motion will be presented and then the oscillator 
model of Hartlen and Currie [87] will be discussed with reference to 


the experimental data presented in the previous chapter. 


B. . EQUATION OF MOTION AND COMPARISON OF RESULTS 
The equation of motion for an elastically-mounted and linearly- 


damped cylinder may be written as 


MX+2Mu ch+ Ko = F(t) (39) 
where X represents the transverse displacement of the cylinder; M, the 
mass of the cylinder plus the displaced mass; c, the damping ratio in 

still water; Wa? the natural frequency of the synchronized cylinder; 

Kos the spring constant; and F(t), the exciting force. F(t) may 

be presented by 


2 


t) = 7edL (C°),.. U (40) 


2 rms =m C (t)/(C,°) 


| rms 


in which Cae is the lift coefficient for the stationary cylinder 
at the corresponding Keulegan-Carpenter number and Reynolds number. 


The ratio C (t)/(C °) 1. represents the amplification of the lift force 
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at or in the neighborhood of perfect synchronization. Even during the 
period of nearly perfect synchronization the exciting force does not 
retain a constant amplitude and frequency. This fact has already been 


demonstrated through the use of the Fourier analysis. In general, it 


4s necessary to consider the harmonics of the lift force. Rewriting 


Eq. (36), one has 
C, (t) = ty C cos(2mnf ,t/T) (36 - repeated) 


Inserting C(t) in Eqs. (39) and (40), one obtains a differential 


equation which may be solved by assuming 
Ge) ay x, cos(2mnf .t/T+9, ) (35 - repeated) 
n= 


The general solution of Eq. (39) is given by 


2 U 
Lake (C0 men 3 
y zl M a rms! ‘SD? Cn! (C are 
a = cs (nwt /THO ) (41) 
: MWe 9 9 2.1/2 a 
eon) ) + (20 ft) } 
n W) 
n 
The relative amplitude X/D is thus given by 
U Ba 
Xx 1 pb D" Bia 5} (Fp) Cf (e Jens 
ies Nae pp nag > 172 (42) 
{{] - ‘oT iit (C2 <G 7 ) 
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Finally, the phase angle reduces to 





-| n 
is cam [ a ] (43) 
jo 
Ww 
n 
The amplification ratios Ch lf C, °) ae have been evaluated for the 


fundamental and two harmonics through the use of the Fourier-decomposed 
exciting force. The results have shown that the said ratio for the 


predominant harmonics at synchronization are given by 


“ip MOL ms SINE yee VS (44a) 
“pg! MOL) ems Soe Oi (44b) 
“pg! OL) ems UR GAG) ee 10) (017 (44c) 
rie ree / (C2) Se jo0 es) 20 (44d) 


in which P, represents the predominant harmonics (see Figs. 114-121). 


The (C, °) values for the corresponding stationary cylinders were 


rms 
taken from the data presented in [101]. 

It is rather remarkable that the amplification ratios at synchroni- 
zation remain fairly constant in spite of the wide range of the test 
Parameters. Consequently, they can be used with confidence in the pre- 
diction of the relative amplitude of the cylinder oscillations. It 


Should be noted, in passing, that there is, at present, no theoretical 


Means to predict the amplitude ratios. The use of the Navier-Stokes 
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equations is limited to very small Reynolds numbers. The discrete-vortex 
model is limited to steady flow about stationary and oscillating cylinders 
[10]. It is not currently possible to apply the discrete-vortex model 

to the hydroelastic oscillation of cylinders in harmonically oscillating 
flow. This is primarily because of the difficulty of the determination 

of the numerous and mobile separation points on the cylinder. Consequently, 
one must rely on the experimentally determined values of the amplifi- 

cation ratios in predicting the cylinder response. 

Equations (44a) through (44c) together with Eq. (41) have been used 
to predict the response of all the cylinders. Figures 167 and 168 show 
two of these plots where the predicted and measured responses are seen 
to compare reasonably well. Similar comparisons with all other cylinders 
proved to be equally satisfactory. 

A critical examination of the relative significance of the second 
and third harmonics on the cylinder response has shown that the response 
is largely determined by the predominant harmonic. This has already been 
anticipated because of the fact that the second and third harmonics of 
the relative displacement are an order of magnitude smaller than the 
predominant harmonic (see Figs. 106-113). Assuming that the fundamental 
harmonic of the exciting force is representable by (Cie and that 


Se = on at synchronization, Eq. (42) reduces to 


Z 
(U,/F 0) (Cams 


16n Mz/{cLD (C,°) 


Ue ) el (45) 


OI ><| 


rms! 
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Fig. 167 Measured and predicted response for a 3 in. smooth cylinder, 
ee 0.\05c. 
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A comparison of Eqs. (44a) and (44d) shows that 


(Cums! (Cems = 6, /(C,°) (46) 


rms LP, rms 


Thus, 1t 1S reasonable to expect that Eq. (45) will yield the fundamental 
harmonic of the relative displacement fairly accurately. As noted 
earlier, higher order harmonics of X/D are an order of magnitude 

smaller and certainly quite negligible. 


Equation (45) shows that a relationship should exist between X/D 


and the response parameter 


M 
Ro = a0. Sma oa (47) 
pv 1 (Cy lane 


for both smooth and rough cylinders. Figure 169 shows the relationship 
between X/D and Ry for all test cylinders. The relationship between the 
two parameters is rather remarkable in view of the fact that only the 


fundamental harmonic of the exciting force has been used and it was 


He 


assumed that We = Ww =w . Figure 169 also shows that the response of 


n nw 
a given body to a given flow must necessarily depend on the starting 
dynamical conditions, i.e., on the dynamics of the same flow past the 
Same body when the latter is held stationary. The use of a response 
parameter given by Mc/(pLD°), as done previously by others [3], without 
er ) ams could not have correlated the response for all smooth and 
rough cylinders. In fact, initially, the experiments were conducted 
with smooth cylinders. It was supposed that roughness will complicate 


the matters further. Contrary to this supposition, roughness and the 
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consideration of the starting dynamical conditions helped to delineate 
the most important parameter governing the hydroelastic oscillations of 
cylinders. 

As noted in Chapter I, researchers have used different definitions 


' for the response parameter. One such definition is 





This particular parameter has been extensively used in correlating 
| the data obtained in steady flow about oscillating cylinders. A 
semi-empirical analysis by Iwan and Blevins [1,89] related the relative 


amplitude in steady flow at perfect synchronization to A. by 


0.44 9,3 + 4.52 $_jV/2 (49) 
A. + 11.94 S A. + 11.94 S 


o|><| 


Griffin, Skop, and Ramberg [61] proposed a least-squares fit to 


the existing data as 
A= raat) + 0.43 a1? (50) 
Finally, Sarpkaya [10, 73] theoretically obtained 
2 ji/2 (51) 


ia 5 211C © 
en OC) in e/ 11-94 + A /KC @) 


in which AUG ea oe is related to Re by 


niZ 





‘, 
SNRs eS) (52) 


° p 
(C) Ven 


Thus, the variation of X/D with R or A/(G °) in harmonic flow 
may be compared with that in steady flow through the use of Eq. (51). 
Figure 170 shows such a comparison for all the data obtained in the 
present investigation. It is rather surprising that not only the 
response in steady flow but also that in harmonic flow is fairly well 
represented by Eq. (51). It is of importance to note that in plotting 


the data shown in Fig. 170 (C,°) values appropriate to each cylinder 


| rms 
| have been used [see Eq. (45)]. 

In summary, it has been shown that the synchronization in harmonic 
flow occurs over a range of U_/ FD values from 5.0 to 7.5 (see Figs. 
Mol, 95, 59, 63), with perfect synchronization at U/ FD ~ 546: The 
amplification factor for the predominant harmonic is about 2.75, and 
the response is a unique function of the response parameter Ro or 


A /(C, °) , for all smooth and rough cylinders. The effect of the 


rm 
Keulegan-Carpenter number and the Reynolds number on the cylinder 
response 1s imbedded in the dependence of CC on the starting 


dynamical conditions, i.e., on the dependence of (C° eae on K and Re. 
Finally, the predominant harmonic of the response is an order of magni- 
tude larger than those of the higher order harmonics, which enable one 

to determine the response through the use of a relatively simple analysis. 
The modulations in amplitude and phase of the response is quite well 
predicted both by a generalized solution based on the Fourier decomposi- 


tion of the exciting force and by a straightforward application of the 


Duhamel 's integral. 
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C. APPLICATION OF THE WAKE-OSCILLATOR MODEL 

several mathematical models have been proposed [86-91]; (see also 
(2, 55, 7/1, 92] for additional models) in an attempt to simulate and/or 
explain some of the experimental results summarized in Chapter I. These 
models do not include the analysis of the flow field and the fluid- 
mechanical] arguments invoked in their evaluation are not altogether 
convincing. Thus their worth should be measured not so much by their 
capacity to obtain functional relations among significant parameters 
that lead to the basic understanding of the phenomenon but by their 
ability to produce results which are qualitatively similar to those 
obtained experimentally. 

The most noteworthy among the oscillator models is the one proposed 
by Hartlen and Currie [87] where a van der Pol-type soft nonlinear 
oscillator for the lift force is coupled to the body motion by a linear- 
dependence on cylinder velocity. This model is based partly on a Sugges- 
tion by Birkhoff and Zarantonello [23], mentioned earlier in connection 
with the estimation of the Strouhal number, and by Bishop and Hassan 
[56] in connection with their now classic experiments with oscillating 
cylinders in uniform flow. The model has its roots in mechanics and 
electricity rather than in the equations of fluid motion. 


The pair of equations which result from this concept are 


ee ° _ 2 
Me ee CG A + x Seeds 2, C (53) 
c -aQ e + ee C 3 + 2 ¢ Cc, = BX (54) 
E oye gl 2 L O L 
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in which 


= dX 7 
X. k/D ’ K. ~ dt ’ T = wit 
DL 
a Se ee eee 
O 9 2.2 
To. 
O 
3 + FD/V (55) 
2, = gli 


where FE is the frequency of vortex shedding by a stationary cylinder 

in steady flow of velocity V and 36 is the corresponding Strouhal number. 
The parameters a and y are the van der Pol coefficients and 8 is the 
interaction parameter. 

In Equation (54), first and fourth terms generate a simple harmonic 
oscillator of normalized frequency QO the second term (the so-called 
negative lift damping) provides the growth of Cs and the third term 
prevents unlimited growth. The fifth term provides the feedback from 
the body motion to the fluid motion and hence to the lift force. 
Equation (54) is not the only form which could provide the desired 
Variation of Cs but it is the simplest. 

Of the three parameters (a, y and B), two must be chosen to provide 
the best fit to the data. This is not too unusual and all other models 
require some retrofitting to experiment. The ratio a/y is related to 


the amplitude of the steady-state oscillation of the lift force by [87] 


C,° = v4q/3y (56) 
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which is the solution of Eq. (54) for 


Experiments in steady flow have shown that [56] in large-amplitude, 
steady-state, vortex-induced oscillation, the displacement and 
exciting force have nearly sinusoidal forms and oscillate at the same 


frequency Wee close to wa? with a phase angle o>. Thus 
X= X sin(Q t) (57) 
and 


eC 


L sin(Q 1 a o) (58) 


Substituting Eqs. (57) and (58) in Eq. (53) and equating to zero the 
resulting coefficients of sin(Q tT) and cos(.t) ; Hartlen and Currie [87] 


obtained 


2 


leas 


> 2 
) x = a, 2, C cos ¢ (59) 


-_ 7 2 : 
ee = a, 2, Cc, sin ¢ (60) 


Likewise, substituting Eqs. (57) and (58) in Eq. (62), dropping terms 


arising from a, and using Eqs. (59) and (60), one has 


Lie. 





2 
{2 ae Se (1 ss 9° i y~l/2 (61) 
0 C ane 2 ee 
(1-2. )* + Gz on 


which shows that 2. <Q orf <f., 
C 0 C O 


Hartlen and Currie applied their analysis to a comparison with 
Jones, et al., [17], data (Re = 3.6 x 10° to 1.9 x 10’, X. = 0.08) 
with good qualitative results. However, some doubt has been raised 
[8] as to whether the oscillations had reached synchronization and 
whether an amplitude of x. = 0.08 was sufficiently greater than the 
threshhold amplitude [44, 62] necessary to bring about synchronization. 
Furthermore, Jones, et al, data did not show an increase in drag in the 
transcritical range of Reynolds numbers. 

Hartlen and Currie also used Ferguson's [63] data for the circular 
cylinder and obtained the results shown in Fig. 171. These results 
are in qualitative agreement with those shown in Fig. 3. There are 
however, some important quantitative differences: 

1. In Feng's data (Fig. 3) the cylinder is seen to continue to 
oscillate at W. = w, on both sides of the region outside the lock-in 
range. This feature is not presented in Fig. 171. 

2. The force and response maximums occur in Fig. 3 at about the same 
V.- Value. In Fig. 171, Cc, remains constant over a broad range of 
Synchronization, reaching a maximum at a V7 value considerably smaller 
(relative to the width of the lock-in range) than that for Xx The 
explanation of these discrepancies and the relationship between the 


particular changes in the wake and/or the elastic system which trigger 


the hysteresis remained unresolved. 
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Fig. 171 Response characteristics in steady flow 
as predicted by Hartlen and Currie model [87] 
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Fig. 172a C. 


by Hartlen and Currie model. 


and x. versus U, in harmonic flow as predicted 
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In spite of its shortcomings, Hartlen and Currie model has been 


| applied to the hydroelastic oscillations of cylinders in harmonic 


flow through the use of Eqs. (53), (54) and 
C(t) =~ “Lp, Sin(Q + o) (62) 


The parameters a and B have been chosen to provide the best fit to the 
data. Representative results for a 4 in. rough cylinder are shown in 
Figs. (172a) and (172b). Apparently, the judicious selection of the 
two independent parameters provides a fairly good agreement between the 
data and the predictions of the model near synchronization. Calculations 
with other cylinders yielded similar results. 

As noted earlier, this model does not include the analysis of the 
flow field and gives no clues regarding the values of a, y, or B. 
Nevertheless, it is worthy of consideration because of its ability to 
predict the general trend of the data once a and B are chosen through 
the use of two data points on the response curve. Additional attempts 
made to improve the oscillator model by making BX term in Eq. (54) a 


non-linear function of ‘a will not be presented here. 
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VI. CONCLUSIONS 


The theoretical and experimental investigation of the hydroelastic 


. oscillations of smooth and rough circular cylinders in harmonic flow 


warranted the following conclusions: 

1. An elastically-mounted cylinder may undergo synchronized oscilla- 
tions when the reduced velocity U. - U/t is in the range of 5.0 to 
a). 


2. Perfect synchronization, at which the response is maximum, occurs 


| wat U. = 5.6. 


3. Synchronized oscillations occur at an average Strouhal number of 
16. 

4. In the region of synchronous oscillations, the rms of the lift 
coefficient and the amplitude of the predominant harmonic of the 
normalized lift force are amplified by a factor of about 2.75 relative 

to that for a stationary cylinder in harmonic flow at the corresponding 
Keulegan-Carpenter number and Reynolds number. 

5. The Fourier and spectral analyses of the exciting force and response, 
in the resonance region, show that the exciting force is well represented 
by its three most important harmonics and the response, by the predominant 
harmonic alone. 

6. The analysis relating the governing parameters predicts fairly 
accurately the response of both the smooth and rough cylinders in the 
Synchronous region. 

7. The measured instantaneous values of the cylinder response are in 
excellent agreement with those predicted through the use of the Duhamel 's 


integral. 
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8. The wake oscillator model of Hartlen and Currie, as applied to 
harmonic flow about transversely-oscillating cylinders, predicts fairly 
accurately the general behavior of the force and response data provided 
that the two free parameters of the model are chosen judiciously. 

9. Outside the synchronization region the cylinders continue to 
oscillate, with very small amplitudes, at nearly their natural 
frequency. 

10. The analyses and the data presented herein are sufficient to 
predict not only the range of occurrence of hydroelastic oscillations 
of elastically-mounted cylinders in harmonic flow but also to evaluate 


the amplitudes of the exciting force and the response. 
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APPENDIX A 


PEAK RELATIVE AMPLITUDES AS A FUNCTION OF K, Un. Fe/ fiw? and are 
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APPENDIX C 


PEAK LIFT COEFFICIENTS AS A FUNCTION OF K, Uns ee and Fe/F. 
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ROOT-MEAN-SQUARE VALUES OF THE LIFT COEFFICIENT 
AS A FUNCTION OF K, Un ere and Fe/f. 
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Fig. 128 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
K = 36 
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Fig. 129 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
K = 39 = 
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Fig. 130 Spectral analysis for a 3 in. smooth cylinder, z = 
= 43 
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Fig. 131 Spectral analysis for a 3 in. smooth cylinder, zc = 0.052, 
K = 46 Meare: 
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Fig. 132 Spectral analysis for a 3 in. smooth cylinder, ¢-= 0.052, 
K = 48 
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Fig. 133 Spectral analysis for a 3 in. smooth cylinder, ¢-=--0.052, 
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Fig. 134 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
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Fig. 135 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
K = 60 
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Fig. 136 Spectral analysis for a 3 in. smooth cylinder, ¢-= 0.052, 
K = 64 
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Fig. 137 Spectral analysis for a 3 in. smooth cylinder, zc = 0.052, 
K = 69 
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Fig. 138 Spectral analysis for a 3 in. smooth cylinder, ¢ 
K = 88 
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Fig. 139 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
K = 99 | 
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Fig. 140 Spectral analysis for a 3 in. smooth cylinder, ¢ = 0.052, 
K = 112 
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Fig. 142 Spectral analysis for a 3 in. smooth cylinder, c = 0.052, 
K = 134 
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SPECTRAL ANALYSIS FOR A 4 in. ROUGH CYLINDER FOR VARIOUS VALUES OF K 
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Fig. 144 Spectral analysis for a 4 in. rough cylinder, ¢ = 0.030, 
K = 29 
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Fig. 145 Spectral analysis for a 4 in. rough cylinder, ¢ = 0.030, 
K = 32 
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Fig. 147 Pk analysis for a 4 in. rough cylinder, c = 0.030, 
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Fig. 149 Spectral analysis for a 4 in. 
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Fig. 152 Spectral analysis for a 4 in. rough cylinder, ¢ = 0.030, 
K = 46 
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Fig. 153 Spectral analysis for a 4 in. rough cylinder, ¢ = 0.030, 
K = 5] 
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Fig. 154 Spectral analysis for a 4 in. rough cylinder, ¢ = 0.030, 
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